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Solutions to equilibrium sequences of irrotational binary polytropic stars in Newtonian gravity 
are expanded in a power oi e = ao / R, where R and ao are the orbital separation of the binary system 
and the radius of each star for R = oo. For each order of e, we should solve ordinary differential 
. . . equations for arbitrary polytropic indices n. We show solutions for polytropic indices n = 0.5, 1, 1.5 

t^r^ I and 2 up to orders. Our semi-analytic solutions can be used to check the validity of numerical 

. solutions. 
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I. INTRODUCTION 



Coalescing binary neutron stars (BNSs) are considered to be one of the most promising sources of gravitational 
> ■ waves for laser interferometers such as TAMA300 §, GEO600 @, VIRGO § and LIGO If we have accurate 

, theoretical templates of inspiraling phase of BNSs, we can determine the mass and the spin of neutron stars from 
the gravitational wave signals in the inspiraling phase We may also extract the informations on the equation of 
^ state of a neutron star from the signals in the premerging phase 0. For this purpose it is important to complete 
theoretical templates of gravitational waves in the premerging phase as well as in the inspiraling phase. Moreover, 
' the study of BNSs in the premerging phase is motivated by the need to provide the realistic initial condition for the 
. merging phase simulations p|-p^. 

In order to obtain accurate theoretical templates of gravitational waves from BNSs around the premerging phase [0 , 
a general relativistic hydrostatic problem with compressible equation of state must be numerically solved. In solving 
this problem, we suppose that BNSs reach quasi-equilibrium states because the timescale of the orbital decay driven 
O ' by the radiation reaction is much longer than the orbital period of BNSs until their innermost stable circular orbit 
(ISCO). (The formalism for solving quasi-equilibrium figures of irrotational binary neutron stars in general relativity 
^ ■ is given in the references p2|-p^.) Furthermore, we regard the internal state of a neutron star as an irrotational or 
• • , nearly irrotational one because the viscosity of a neutron star is not large enough for synchronization even near the 
ISCO 

Recently Bonazzola, Gourgoulhon and Marck have numerically studied the irrotational binary neutron stars in 
5J] ' general relativity with the conformally flat condition [p^|-^ . More recently, Uryii and Eriguchi have also numerically 
[ solved the same problem as that of Bonazzola, Gourgoulhon and Marck . In order to check the validity of these 
numerical calculations, it is necessary to compare them with analytic or semi-analytic ones. However, we have not 
obtained such analytic solutions in general relativity yet. There are solutions of irrotational BNSs only in the first 
post-Newtonian (IPN) approximation of general relativity. For example, Lombard!, Rasio and Shapiro have semi- 
analytically studied irrotational BNSs by using the energy variational method |2^, and one of the authors of the 
present paper (KT) has analytically calculated equilibrium sequences of irrotational BNSs by using the tensor virial 
method [ p3| . The former one restricts the internal motion within the plane orthogonal to the rotational axis assuming 
the shape of the star to an ellipsoidal one in order to treat compressible equation of state. While in the latter one 
the velocity component along the orbital axis is included at IPN order and the shape of the star at IPN order is not 
restricted to an ellipsoidal one although the fluid is incompressible. 

The situation is also the same in Newtonian gravity, that is, even in Newtonian gravity useful analytic or semi- 
analytic solutions for equilibrium sequences of irrotational binary systems are not obtained. Let us consider numer- 
ically constructed stationary structures of irrotational binary stars computed by Uryii and Eriguchi in Newtonian 
gravity j2^j2^. One may think that semi-analytic solutions by Lai, Rasio and Shapiro ||2^] may be used to check 
the validity of numerical solutions. However, in numerical solutions of Uryii and Eriguchi the velocity component 
along the orbital axis exists while in those of Lai, Rasio and Shapiro such a component is assumed to be zero from 
the beginning. Therefore new analytic or semi-analytic solutions are needed to check numerical solutions even in 
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Newtonian gravity. Such a check of numerical solutions is extremely important because in the numerical calculation, 
there is a possibility to obtain another solution although the binding energy of a binary neutron star is almost the 
same value, and to lead a different conclusion 

In the previous paper |29| , we showed such a new, almost analytic solution to an equilibrium of irrotational binary 
polytropic stars for the polytropic index rt = 1 in Newtonian gravity by expanding all physical quantities in a power 
of e = aa/ R, where R and oq are the orbital separation of the binary system and the radius of each star for R = oo. In 
that paper, we have extended the method developed by Chandrasekhar more than 65 years ago for corotating fluids 
pO| , pl| to the one for irrotational fluids. In this paper, we show semi-analytic solutions for arbitrary polytropic indices 
by numerically solving ordinary differential equations. 

This paper is organized as follows. In §2, we formulate the method to solve the irrotational binary polytropic stars. 
In §3, the physical values, i.e., the central density of a star, the orbital angular velocity, the total energy and total 
angular momentum of the binary system are calculated. The numerical results are presented in §4. Section 5 is 
devoted to summary and discussions. 

Although a binary system consists of two stars, we pay particular attention to one of two stars. We call it star 1 
whose mass is Mi and the companion one star 2 whose mass is M2. In this paper, we adopt two corotating coordinate 
systems. First one is a Cartesian coordinate system X whose origin is located at the center of mass of the binary 
system. For calculational convenience, we choose the orbital axis as X3, and we take the direction of Xi from the 
center of mass of star 2 to that of star 1. The second coordinate system is the spherical one r = (r, 9, ip) whose origin 
is located at the center of mass of star 1. We use units of G = 1. 



II. FORMULATION 



Since we treat irrotational fluids in Newtonian gravity, the basic equations for obtaining equilibrium configurations 
of binary systems are the equation of state, the Euler equation, the equation of continuity and the Poisson equation: 



P = Kp^ 



dP 



-U + -v^ - t) • (O X r) 

V • V = -(v - Clxr) ■ — , 
P 

AU = -47rp, 



= 0, 



(2.1) 
(2.2) 

(2.3) 
(2.4) 



where P, p, n, U and Q are the pressure, the density, the polytropic index, the gravitational potential and the orbital 
angular velocity, respectively. K is a constant related to entropy and v represents the velocity field in the inertial 
frame. The gravitational potential U is separated into two parts, i.e., the contribution from star 1 to itself U^^^ and 
that from star 2 to star 1 U^^^. We can also express the gravitational potential contributed from star 1 to star 2 as 
jji-*2 ^ These gravitational potentials are written as 
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where the superscript ' means the term concerned with star 2 and Pi denotes the Legendre function of order I. Pn 
is the reduced quadrupolc moment defined by 
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In the irrotational fluid case, we can express v as a gradient of a scalar function i.e., 



V = V$. 

Following the Lane-Emden equation we first express the density as 
with ^ — r I a. Here the definitions of o. and a' are 



a i 



1/2 



An 



1/2 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



with pc being the central density. Let us define e = ag/R, where R and oq are the orbital separation of the binary 
system and the radius of star 1 for R = oo The radius of a spherical star gq is given by oq = a£,i as usual. 

Also we define e' = a'^/R for star 2. Note that in the case of Newtonian gravity, we have two degrees of freedom for 
describing the differences of mass, i.e., pc and a. Therefore, we can express the differences of mass by changing only 
Pc with fixing a. In this case, we can make the radius of each star and each expansion parameter coincide. Although 
we have adopted such a situation in Ref. p9| ], in the present paper we express equations concerned with star 2 by 
using e', a' and so on. 

We expand in a power series of a parameter e a^ 



e 



oo 

E 

i=0 



(2.12) 



Since the shape of star 1 is spherical when R is large, the lowest order term of Q is the solution of the Lane-Emden 
equation. Then we expand O.^ by spherical harmonics as 



Now we consider the orbital motion of star 1. In the spherical coordinate system, it becomes 

nxr = nR{n x ^)^^,^ + Qa{n x ^)^^^, 

where 

1 



1+P 

(0, 0, ^sin 



(sin 6* sin (p, cos 6* sin (p, cost^). 



(2.13) 

(2.14) 

(2.15) 
(2.16) 



and p = M1/M2. The first term on the right-hand side of Eq. ( p.l4[) comes from the orbital motion of the center 
of mass of star 1 and the second term comes from the fluid motion around the center of mass of star 1. The orbital 
motion of star 2 becomes 



nxr'^ nR{n x + na'{n x ^)^^_^, 



where 



X €)lr6 = ^Y^^^^'^^'^"^'^'' cosei'sin^', costp'), 
{nx^y^^^ = iO, 0, ^'sin^'), 



(2.17) 

(2.18) 
(2.19) 



^Note that we express the Lane-Emden function for star 2 as and expand it in a power series of a parameter e' as 



3 



in the spherical coordinate system whose origin is located at the center of mass of star 2. Here we have expressed the 
spherical coordinate system for star 2 as r' = ( r' , 6 ' , cp') with ^' = r' /a' . 
Next, we rewrite the equation of continuity (^.3|) as 



A$ = -^(V* - r2 X r) 
The condition for <I> at the stellar surface is 

(v$ - rj X r) • (ve)[ 

since 6 = at the surface. We expand $ also as 



ve 
e ■ 



surf 



(2.20) 



(2.21) 



(2.22) 



i=0 



When R is large, the shape of star 1 is spherical and star 1 has only the orbital motion of the center of mass in the 
inertial frame with no intrinsic spin so that the gradient of the lowest order term of $ should agree with the orbital 
motion of the center of mass of star 1, i.e.. 



V<I>o = QRin X $,) 



orb ■ 



Considering we normalize $ as 



We again expand $i by spherical harmonics as 



where Y/^ is defined by using the spherical harmonics Y/^ as 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



with constants Ci„i- Note that the gradient of the lowe st order term of $ has order except for the implicit 
dependence in fi. When we normalize $ using Eq. (2.24), the lowest order term of $ has order e°. We must keep 
this difference in mind. 

The orbital angular velocity is derived from the force balance equation, i.e., the first tensor virial relation, defined 
by M 



star 1 



dP 
dxi 



where xi — rsinO cos ip. This equation is expressed by using Eq. (2.2) as 



= 



d-^xp 



dxi 



(2.27) 



(2.28) 



If we substitute Eq. (2.14), 8o and <f>o into Eq. (2.28), we obtain the orbital angular velocity in the lowest order as 

Mtot 



^2 - 
^^0- ^3 



where Mtot — Mi + M2. Note here that we also expand as 

00 



(2.29) 



(2.30) 
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Finally, we express the gravitational potential by rewriting Eq. (2.2) as 

U = K{1 + n)p^ + - w . (n X r) + C/q, 



(2.31) 



where [/q is constant. Substituting Eq. (2.31) into the Poisson equation (2.4), we obtain the equation to determine 
the equilibrium figure as 



a^Ae = -9" 



1 



(V$)2 - 2(V$) • ($7 X r) 



(2.32) 



Now a solution can be obtained iteratively. Firstly, Oj is determined by demanding that the gravitational potential 
and its normal derivative are continuous at the stellar surface |30|, that is. 



C/™t(e = S) = C/e,t(e = S), 



where ^(6*, (p) expresses the surface, 



e(e = s(0,^)) = 0. 



(2.33) 
(2.34) 



(2.35) 

2 



Substitut ing ^ and Eq. (2.14) into Eqs. ( 2.2Cl| ) and ( 2.2S| ), we obtain $i and ilf . Secondly we substitute <i>i and Q 
into Eq. (2.32) and derive Qi+i- We have continued this procedure up to order in this paper. 



III. DERIVATION OF THE PHYSICAL VALUES 

In this section, we calculate physical values such as the baryon mass of star 1, the orbital angular velocity, the total 
energy and total angular momentum of the binary system. We will just show the main results. All the details will be 
shown in the Appendix. 



A. Mass and Central Density 

The baryon mass of star 1 is calculated by integrating the density over the volume of star 1 as 
Ml = / d^xp, 

J star 1 



= / d^d9d<pf sine Q'^^ + e^nQ'^-^Qs + e-^n&^-^Oi + e^n&^-^Qs 

4e^'06 + ^(n-l)e^^ei}]. 



-e^n< 



Note here that the integration range for ^ is [0, ^{9, (p)]. 
For n > case, the mass of star 1 becomes 



1 



77TZr-7^'^®o (Ci) + ^Vo,c(Ci 



where we define 



.eo,c(a)^V 1080,^(6) 
- -4-Pca^^?eo,(er) [l + ^^(^^e^) (ca + ^^^0(6)) 



(3.1) 



(3.2) 



(3.3) 
(3.4) 

(3.5) 



We have used Eq. (|H43| ) in deriving Eq. (3.4) from Eq. (3.3). Here cg is a constant, and one can see the expression 
of ce in Appendix ^We assume that the baryon mass Mi — — 47rpcoQ!oCi0o,c (^i) is conserved throughout the 
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equilibrium sequences of the binary system, where pco and ao denote the values for a spherical star. Therefore, the 
central density and a are expressed as 



where we define 



Pc = pco ( 1 - e 



6 ^Pc 



K{1 + n)p-o 



2n pc 



Pc 3 - n\$,iQo,d^i)^ ^ ' 



(3.6) 
(3.7) 

(3.8) 



B. Orbital angular velocity 

We obtain the orbital angular velocity from Eq. ( ^.28 ) as 



Mu 



2i?2 

9e2 



/3ln 
V Ml 



O(e^) 



ag L 2ag V Ml M2 / ^ ^ 



where J^n denotes ^fn/e^. 



The total energy is written as 



C. Total energy 



E = litot + {Wself)tot + {W,nt)tot + TtoU 



(3.9) 



(3.10) 



where litot, {Wseif)tot, {Wint)tot and Ttot denote the total internal energy, the total self-gravity energy, the total 
interaction energy and the total kinetic energy, respectively. We calculate each energy in the following. 



1. Internal energy 

The definition of the internal energy of star 1 is 

ni=n f d^xP. 

For 71 > case, we obtain 
Ml 



(3.11) 



Hi = 



l + naoe?(eo,c(6))' 



1 - e' 



; 5 — 71 Spc\ 



2". Pc I Jo 



(3.12 



Therefore, the total internal energy is expressed as 



^The expression of Spc/pc seems to diverge if n < 1 when one see Eq. (3.3). However physically Spc/pc should be finite. In 
reality we determined Spc/pc by using the virial relation as shown in the Appendix^ 
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+ 



n 



Ml 



1 + n' a'^^[-\eo.,^ic[)y 



(3.13) 

where the superscript ' means the terms concerned with star 2, and ~ also means the functions concerned with star 2. 



2. Self-gravity energy 



The definition of the self-gravity energy of star 1 is 

1 



w;. 



For n > case, we obtain 
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where fi is defined as 



M = 3^ = -( )60o,«(6)- 



The self-gravity energy of star 2 is 

1 + n' 
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Here, Cg is a constant similar to ce, and 



Therefore, the total self-gravity energy becomes 
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(3.14) 
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3. Interaction energy 
The definition of the interaction energy of star 1 is 

Wint,l = -\ I (t'xpU^^\ 
^ J star 1 



M1M2 _ 33M1M2/3J21 ,3fn\ 



2ao 



Therefore, the total interaction energy becomes 



, M1M2 33M1M2 



ao 



2ag r Ml M2A 



(3.20) 
(3.21) 

(3.22) 
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4- Kinetic energy 



The definition of the kinetic energy of star 1 is 



^ J star 1 
M1M2 



2ao(l +p) 

The kinetic energy of star 2 is also given as 



9f_ 

2al 



\ Ml 



^11 , g/3^JJ_ 

Ah 



To = 



M1M2P 
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Then the total kinetic energy becomes 
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2ao 



M1M2 f 3 ^11 ,3 ^11 
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(3.23) 
(3.24) 

(3.25) 

(3.26) 



D. Total angular momentum 

The total angular momentum is calculated from the equation; 



J 



d^xp{R X v), 



star 1, star 2 



where R for star 1 is written as 



and for star 2 is 



R 



R' = 



R 



1+p 



(sin 61 cos cos 61 cos — sincp)+r(l, 0, 0), 



pR 
1+p 



(sin 6*' cos (p', cos 6*' cos (/j', — sm(p') + r' {1, 0, 0) 



(3.27) 



(3.28) 



(3.29) 



in the spherical coordinate system. Note that since the velocity field v — V$ is written in the spherical coordinate 
system, it is convenient to calculate R x v in the spherical one. However, we should be careful to integrate it over 
the volume of each star, because in the spherical coordinate system, the unit vectors depend on the coordinates. For 
example, the unit vector (1,0,0) in the Cartesian coordinate system becomes (sin0 cos cos cos — sin (/s) in the 
spherical one. 

As a result, we obtain the total angular momentum; 



J = (0, 0, J) 

in the Cartesian coordinate system. Here J is calculated as 

M1M2 



J 



Ml + M2 



-R'^n 



1 + higher order terms than 0(e 



(3.30) 



(3.31) 



E. Virial Relation 

By using the above energies, we can calculate the virial relation which is written as 

-Hi + ^H2 + {Wself)tot + {W^nt)tot + 2Ttot = 0. (3.32) 

n n' 

This equation is used to check the solutions as shown in the previous letter [ p9[ . However in this paper we use it to 
determine Spc/pc (see Appendix |). 
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IV. NUMERICAL RESULTS 



In the following sections, we set Mi — M2 — M, n — n', a — a', pc — p'c and so on, i.e., the identical star binary 
for simplicity, and investigate the cases of different polytropic indices. 

First, we show the Lane-Emden function of star 1 up to O(e^). It is written as 



6^94 + £'^65 + £^66, 



where 



63 

94 

65 



(3V2(0-P2(sin6'cos(p), 

('*V3(C)'P3(sin6'cos(p), 

(^V4(0-P4(sin6'coscp), 

^^Vo(f) + '^^V2(C)-P2(sin0cos 

+(^V4(^)^4(sin6'cosc/?) + (^V5(0^5(sin6lcosv3). 



The functions Oq, and their derivatives for each polytropic index are given in Tables H/ml 
Next, we give the velocity potential up to O(e^). We can write it as 



$ = 
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(4.1) 



(4.2) 
(4.3) 
(4.4) 

(4.5) 



(4.6) 

(4.7) 
(4.8) 
(4.9) 

(4.10) 



The functions ^'*V2, ^^^3 and ^^■'04 and their derivatives are shown in Tables IX -XII. The derivatives are related 
to the velocity field in the inertial frame. Note here that in Eq. (4.7) we write the dependence of the mass ratio 
p = M1/M2 although we consider only p — 1 case in this section. 

Here we directly show the existence of the velocity component along the orbital axis. In this paper, we take as 
the orbital axis. Since the origin of the corotating coordinate system X is located at the center of mass of the binary 
system, it is convenient to consider another Cartesian coordinate system which origin is located at the center of mass 
of star 1. We call it a; = {xi, X2, X3). The relations between x and the spherical coordinate system r are 



Xi 
X3 



r Sill tl COS (p, 
r sin sin 93, 
r cos 6, 



as usual. Then, the velocity component along the orbital axis in the inertial frame is written as 

■i^3 



8x3 ' 



In the corotating frame, we can express the velocity field u as 



-na{n X I) 



fig 



e3V$4 + e V$5 + e V$e 



(4.11) 
(4.12) 
(4.13) 



(4.14) 



(4.15) 
(4.16) 
(4.17) 



where V = aV is a nabla defined by using {£^,9,ip) and given in Eq. (A4). Then, the velocity component along the 
orbital axis in the corotating frame becomes 
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M3 — riflo 



where we have used the normahzed coordinate system, 

x\ — ^ sin Q cos 
X2 = ^ sin sin 

^3 = C COS0, 



(4.18) 



(4.19) 
(4.20) 
(4.21) 



for the calculational convenience. Note that since V$ is proportional to fi/e, the equation ( 4.16 ) does not express 
the real dependence on e. 



Next, we calculate the term for $4 in Eq. ( 4.18 ). Then we obtain 

a$4 _ „9$4 sin6'a$4 
^ -cos^^-— 

= 3 sm 6* cos 6* sm 2(p. 



(4.22) 



If *^'*V2(0 is proportional to the velocity component along the orbital axis disappears^. However, '^^■'(/)2(C) is not 
proportional to ^ at all for n > 0. Moreover, we can show the velocity components along the orbital axis for $5 and 
$6 as 



dx-i 



= 6 
= 30 



id^%^ 4(6)04 



, , (1 - 5sin^ 6'cos^ (^) , 



(3 — 7 sin Q cos (p) 



6(6)04 



smf/cost'smi^, 



sm t'cos6'sm</?cos(^. 



(4.23) 
(4.24) 



In these equations, we can see that even if (^Va oc ^ and (6)04 cx: there remain the velocity components along 
the orbital axis. The velocity component along the orbital axis is order higher. However, if one would like to pay 
attention to the internal state of neutron stars in the binary system, it should be taken into account. 
We give also the other components of the velocity field as 



,3^ 

dx\ 



5 9<f6 



dx\ dxi 



U2 



^3 9$4 _j_^4^ 

dx2 dx2 



dx2 



where 

a$4 

dxi 



= sm cos ip- 



9<I>4 



COS COS (f d^4 sin ip 9$4 



= 6 



/d(4)02 2(4)02 X . 

[—£ H '''' ' 
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sin^siniyS, 



2 ^ ■ V2 

COS (y5 + 



9$4 
9X2 



sm f sm Lp- 



d^ i J e 

9<I'4 COS 9 sin </? 9$4 cos ip (9$4 





= 6 
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2(4)02 










= 6 


/d(5)03 


3(5^3 


92)1 


V 






= 6 


/d(5)03 


3(5^3 


9^2 







. 2 , 



smWcosiyj, 



(1 - 5sin^ 0cos^ ip) 



sm y sm p cos 



(1 - 5sin^ eicos^ ip) sin2 6'sin^ ip + + sin^ 61(2 sin^ - 5 cos^ (p)| 



(4.25) 
(4.26) 

(4.27) 
(4.28) 
(4.29) 
(4.30) 
(4.31) 
(4.32) 



^For the n = case, i.e., the incompressible fluid case, ('*)02 is proportional to 
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d^6 

dxi 

dx2 



= 30 



30 



^ (3 — 7 sin^ 6 cos^ (p) sin'^ 9 sin cos^ H — 5 sin^ 6* cos^ (^) sin 9 sin 



(3 — 7 sin^ 6* cos^ (p) sm^ 9 siv? p cos ip 



H — ^ 1 3 + sin^ 9{7 cos^ <^ - 6 sin^ v?) | sin 6* cos Lp 



(4.33) 



(4.34) 



In Table XII], we show the first zero point of 0o, the orbital separation at the contact point, the energies, the 



quadrupole moment, and the change in the central density. In this table, we have used 

E 



E = 



E. 



self 



eE, 



'point 



e'E 



quad ■ 



Each energy is given in Table KIV, where we have defined as 



tot — ^^self ^ ^ ^^quad^ 
6/ 



{Wself)tot = (yVself)self + ( W^sei/ )gMod , 
{Wint)tot = ^{Wint)point + {Wint) quadi 
^tot — ^-^point ~1~ € -^o' 



' quad- 



Accordingly, we can express 



Eself — n^ei/ + {Wself)self, 
Epoint — iy^int^point ~^^point: 

Equad — Ilguad ~t~ i^^self^quad {y^int^quad ^quad- 

The total angular momentum for the identical star binary becomes 

M o 

J = —R^n, 



and the orbital angular velocity is calculated as 



1 + 



M 



(4.35) 



(4.36) 
(4.37) 
(4.38) 
(4.39) 



(4.40) 
(4.41) 
(4.42) 



(4.43) 



(4.44) 



We show the total energy, the t otal angular momentum and the orbital angular velocity along the equilibrium sequences 
of the binary system in Table XV. 

In Figs. ^]-^, we show (a) the total energy and (b) the total angular momentum as functions of the orbital separation, 
and (c) the orbital angular velocity as functions of the total angular momentum for each polytropic index, where we 
have used 



Po 



(4.45) 



In these figures, we compare our results (solid lines) with those of numerical calculations computed by Uryii and 
Eriguchi (open triangles) and those of semi-analytic ones proposed by Lai, Rasio and Shapiro (filled circles). 
For the smaller orbital separation, we must expand the physical values up to higher order than O(e^) in order to 
include the effect of the spin of each star. In the case of the irrotational binary system, there is no intrinsic spin for 
the large distant stars. However, when the stars reach close range, the small spins are produced by the deformations 
of the stars. For example, the effect of these spins appears in the total energy at 0(e^)^. The effect of the spin in 
the case of the smaller polytropic index is larger than that in the case of the larger polytropic index. The reason is 



*This is because the spin kinetic energy is given by the volume integral of (e'*V$4)^ 
definition of V is in Eq. ( A.4 ) . 



n'^e'^{V^4f, where Q.'^ is 0{e^). The 
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that the matter i s not concentrated in the case of the smaller polytropic index, and the quadrupole moment becomes 



larger (see Table Kill ). Therefore, at the smaller orbital separation, the deviation between our results and those of 
Uryii and Eriguchi for n = 0.5 (Fig. |^) is larger than that for n — 1.5 (Fig. ||). 

Finally, we represent the accuracy of our numerical calculations by comparing numerical and analytic solutions for 
n — 1 cases pyl. In Fig. a, we show the absolute values of the relative error for (3)^2, (4)^3^ (5)^^^ (6)^2, (6^4 and 
'•^Vsj where we define the relative error by 



Numerical solution 



Analytic solution 



(4.46) 



We can see from Fig. || that the relative errors are within 10 Therefore, we regard that our results have enough 
accuracy. 

V. SUMMARY AND DISCUSSIONS 
A. Summary 

In this paper, we have calculated the equilibrium solutions of irrotational binary polytropic stars in Newtonian 
gravity by expanding all physical quantities in a power of e. We have presented the results of the cases of several 
polytropic indices (n = 0.5, 1, 1.5, and 2). In particular, we have shown the velocity fields by solving the equation of 



continuity. It is found that there exists the small velocity component along the orbital axis (see § [V ) . It agrees with 
the numerical calculations performed by Uryu and Eriguchi [|2^ and Bonazzola, Gourgoulhon and Marck p^jl9| , |3^ . 

Furthermore, we have given the figures and tables of the total energy, total angular momentum and orbital angular 
velocity along the equilibrium sequences for each polytropic indices. We can see from these figures that our results 
agree with those of Lai, Rasio and Shapiro and Uryii and Eriguchi |^ for R/uq > 3. 

Since our solutions are correct if e ^ 1, they can be used to check the validity of numerical solutions. For any 
numerical codes, one can ask to solve an equilibrium for large R, and compare numerically derived velocity distribution 
and so on with our semi-analytic solutions. 

However, since we expanded physical quantities up to O(e^), it may not be enough to discuss about the behavior 
of solutions for small R. In order to apply our solutions in the case of small R and to check the validity of numerical 
codes in this case, further higher order calculations are needed. 

B. Discussions 

It is important to compare the velocity field which we obtain by solving the continuity equation with that given 
by Lai, Rasio and Shapiro. We can see the velocity field they give in their papers p^ , ^ or in the Chandrasekhar's 
textbook I^J]. In the irrotational case, it becomes 

{uLRs)i = 2 / 2 ^^2, (5.1) 
+ 02 

(wLfls)2 = 2 ,^ 2 ^^!' (5-2) 

af + 02 

{uLRsh = 0, (5.3) 

in the corotating frame. Here denotes the length of the principal axis parallel to the x^-axis. We can rewrite these 
component of the velocity field as 

2 2 

{ULRS)1 = ^X2 + (^^-^)nX2, (5.4) 

2 2 

(uLRs)2 = -^xi + ( "^ ~°^ ) nxi. (5.5) 

\ Oil n~ 0"-) ' 



f 

^3 



The second terms in the above equations are order e because the deviation between a\ and 02 is produced by the 
tidal force and the effect of the tidal force is order e^. 
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On the other hand, if we restrict the form of the function in our velocity field as '^''^2 = (&2/Ci)^^> ''^Vs = (Qifcs/Ci)^^ 
and (6)^4 = {aHi/ii}^ as in the case of the incompressible fluid, we can express our velocity field by using Eqs. 
(|^, ( |42^ ) and ( mI ) as 

{Upresent)l = ^^2:2 + ^[Qe'b2X2 - A&t^b'iXiX2 + ^Ot^biX2{r^ ~ Sx^)], (5.6) 

{upresent)2 = ~^xi + Vl [&e'b2Xi + &e%:i{r^ - bxl + 2x1) + 30e^54r(3r2 + Ixl - 6x^)] , (5.7) 
[upresent)z = [l2e'*fe3a;22:3 + 180e'^642:ia:2a;3] , (5.8) 

where = x\+X2+x\. Therefore, we find that the form of the velocity field we obtain in the case of the incompressible 
fluid coincides with that given by Lai, Rasio and Shapiro up to 0{e^) including the value of 62- This means that 
the velocity field of Lai, Rasio and Shapiro is correct only in the case of the "incompressible" equation of state and 
"ellipsoidal" figures, because the velocity field at order is produced by the ellipsoidal deformation of star 1 (see 
Appendix ^). 

Finally, we discuss about the configuration of each star. When we pay attention to star 1, the equation for the 
stellar surface is written as 



(5.9) 



3i!V2(?l) . „ ^ , 4 ^'V3(gl) p , . ^ , , 5 (^^4(^1) o , . ^ , 

e — P2(sin 9 cos ip) + e — - — — — F3(sm0cos(y5) + e — - — -— P4(smt^cos(y5) 



|0o,«(6)l |eo,e(6) 



|0o.e(ei)l 

16^ i^^i^^ + ^(^^)) l35^^(^^"^^°^^) + ^P2(sm.cos^) + -} 
+(^Vo(Ci) + (^^2 (6)^2 (sine cos (^) + (6)^22 (eO^Kcos^) cos 2^ 
+^'^^M^i)Piisin0cosip) + (6)i/,5(^i)F5(sin6icos(^)l , 



(5.10) 



where 5*^ are defined in Appendix ^ and we have used the relation ( |C27 ). Accordingly, we can express the length of 
the principal axis. Although the real length of the axis is written as a2, we show the results divided by a. 

(1) The opposite direction to star 2: 

('^2(6) 



4i!MiI , 5i!MiI 

|eo.c(Ci)l ' ^ |eo,c(a)l |eo,«(6)l 



('¥2(6) /^'V2(ei) , d(^^2 



|eo.«(Ci)|L|0o,c(a)|V 6 
+(6^4(6) + (^¥5(6) 



(a) + ('¥0(6) + ^'¥2(6) + 3(6V22(6) 



(5.11) 



(2) The direction to star 2: 



^^¥2(^1) 4 ''V3(a) , 

|eo,«(a)l ' |eo,«(a)l 
1 r (^V2(ei) / ^^V2(a 



|0o,«(a)i 

^(3^2 , 



|0o,e(a)|L|0o,«(a)l V a 

+(6V4(a)-(^V5(ei 



(6) + ('Vo(a) + ^'¥2(6) + 3(6^22(6) 



(5.12) 



(3) The (opposite) direction to the orbital motion 

"^U' 2/ " "V2'T/' 

3 ^'¥2(6) 



= a - 



3(5)^4(^1) 



2|eo,4(a)l 8190,4(6)1 

1 r ('¥2(6) /('V2(a) 



eo,f(ei)|V" 



|eo,«(a)l L4|eo,e(ei)l V e 



(Cl) 



^'¥0(^1)- 2(*'V2(ei)- 3(^2(6) 



(5.13) 
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(4) The direction parallel to the rotational axis: 



•(0,0) = S(-7r,0), 



+ e- 



3(^V4(ei) 



2|eo,c(a)l 8|eo,«(a)l 
1 r (=^¥2(6) /('V2(a) 



©o,«(a)l L4|eo,e(6)l 



(5.14) 



We can see from these equations and Tables | - VIII that the axis to star 2 is the longest, and the deviation between 
the axis to star 2 and that opposite to star 2 appears at order £■*. On the contrary, the deviation between the axis to 
the orbital motion and that parallel to the rotational axis appears at order e^, and the difference is the effect of the 
deformation induced by the spin of the figure (^^¥22)- 

When we see the quadrupole moments in Eq. (5.10), we find that the coeffient of the higher order term is not large 

as 



|0o,«(ei)l 



M t 

l7|eo. 



(6)1 



) (3V2(a)J 



6 



P2{s\n 9 cos ip). 



(5.15) 



Therefore we can expect the convergence of these terms. However, there is another quadrupole moment induced by the 
spin of the figure. The term e^((^¥22(Ci)/|0o,4 (Ci) 1)^*2 (cos 0) cos2ip seems to be as effective as the leading quadrupole 
term for the smaller orbital separation. This means that the terms concerned with the spin of the figure which appear 
at order in the total energy may change the behavior of the total energy. Furthermore, the hexadecapole moments 
inEq. (pO|), 



|0o,«(a)i 
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35(5^4(6) I0o,aa)iv a 



A.) 



'^¥4(6 ) 

(5V4(C 



^1 



P4(sin6'cos(/3), 



(5.16) 



does not seem to converge at order e^. However, since these terms will appear at order in the total energy, they do 
not have so much effect. Anyway, if we discuss about the behavior of the total energy, the total angular momentum 
and so on for the smaller orbital separation such as R/a^ < 3, we must calculate at least up to order e^. 

ACKNOWLEDGMENTS 

We would like to thank K. loka and K. Nakao for useful discussions. This work was partly supported by a 
Grant-in- Aid for Scientific Research Fellowship (No. 9402: KT) and Grant-in- Aid of Scientific Research (No. 11640274, 
09NP0801: TN) of the Japanese Ministry of Education, Science, Sports and Culture. 

APPENDIX A: SUMMARY OF THE EQUATIONS AND THEIR BOUNDARY CONDITIONS 



In the following appendices, we derive equations order by order of e, and describe in detail. In this appendix, we 
summarize the equations which should be solved numerically with their boundary conditions in each order. 

First of all, we give the total equations which include all terms up to 0(e®). The equations for determination of 
the velocity potential and stellar configuration are written as 



where 



Ae = -9" 



Vl- - (f2 X i)orb - 7-(f^ X Clf.g 



ve 



(V$)2 - 2(V$) • {(fj X £,)orb + ji^ X £,)f,g'^ 



A = a^A 



1 d „ 



92 



e sin^ 9 9(^2 ' 



dC ^99 ^sm9dip 



(Al) 
(A2) 

(A3) 
(A4) 
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The boundary condition for the velocity field is 



v$ - (f2 X ^)ort - -(o X ^) ■ (ve) 

t.1 J 



surf 



= 0. 



The internal and external gravitational potentials which should be matched at the stellar surface are 



-^.5](-l)'.'(f ) V,(sin0cos^) - -^(4i)63[l - 34Fi(sin0cos^) 



/=0 

Z,m /,m Z,m 



where Kq S'Hd Ki:im are multipole moments. Uq is constant and expanded as 

Ua = ca + eci + e^C2 + e^cg + e^c^ + e^cg + e^ce- 
The velocity potential and configuration hmction are expanded up to 0(e®) as 

$ = $0 + e*i + £^^"2 + e^$3 + e'^*4 + £^^$5 + e^^6, 
e = Go + eGi + e^Ga + e^Gg + 6^64 + £^65 + e^eg. 

The boundary conditions for Q are 

where S denotes the first zero point of the function 8, i.e., the stellar surface and is formally expressed as 

E{9, ^) - a + eSi{e, (^) + e'S2{e, (p) + e^S^iO, ^) + e^S^iO, + e>S^{e, ^) + e^'Seie, (^). 
The regularity conditions at the center of each star is 

99, 



9C 



- 0) - 0, 



and also we normalize the central value of Q as 



e(C = o) = i. 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 
(AlO) 



(All) 
(A12) 

(A13) 

(A14) 

(A15) 



In the following, wc show the equation for the veloc ity pot entia l, its boundary condition, the equation for determi- 
nation of the figure and its boundary conditions (All) and (A12). 



The equation for the velocity potential is 



AI>o = —n 



1. 0th order 



Vl>o -{Clx ^)o„ 



60 ' 



(A16) 



and its boundary condition is 
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Vh-i^xOorb -(Veo) 



= 0. 



The equation for determination of the figure is 



1 d / 2ti6o 



and its boundary conditions at the stellar surface are 

Co = 1^, 



— - -72- 



The regularity condition at the center of the star is 
and the normalization of Qq at the center of the star is 

eo(^ = 0) = 1. 



2. 1st order 



The equation for the velocity potential is 



A$i = -n 



V$l - X ^)fig 

?1 



V9o 

eo ' 



and its boundary condition is 

(v$i)-(veo)|j^=o. 

The equation for determination of the figure is 

and its boundary conditions at the stellar surface are 

H{3 + 2p) o , 



Cl 



2(1 ^1 t;r^^' 



The regularity condition at the center of the star is 



= 0) = 0, 



and the normalization of 0i at the center of the star is 

ei(^ = o) = o, 

because we take ©o(^ = 0) = 1. 
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3. 2nd order 

The equation for the velocity potential is 

A$2 = -n(Vi>2) • (A30) 

and its boundary condition is 

(V$2)-(Veo)|^^=0. (A31) 
The equation for determination of the figure is 

AOs = -ne^-^Oa, (A32) 
and its boundary conditions at the stellar surface are 

C2 = -^S, + Y.^^r, (A33) 



The regularity condition at the center of the star is 

992 



= 0) = 0, (A35) 



and the normalization of ©2 at the center of the star is 

62 = 0) = 0, (A36) 

because we take ©o(^ = 0) = 1. 

4. 3rd order 

The equation for the velocity potential is 

A$3 = -n(Vi3) • (A37) 

and its boundary condition is 

(V$3)-(Veo)|^^=0. (A38) 

The equation for determination of the figure is 

AOs = -ne^^Os, (A39) 

and its boundary conditions at the stellar surface are 

C3 - -^P2(sinecos^) = -^S, + Y: "-^rr, (A40) 
1+p ?i , £ 

^ l,m 



The regularity condition at the center of the star is 

d&3 



= 0) = 0, (A42) 



and the normalization of ©3 at the center of the star is 

63 = 0) = 0, (A43) 

because we take Oo(^ = 0) = 1. 
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5. 4th order 



The equation for the velocity potential is 



(vl>4)-(veo)--(i^x$)/,,.(ve3) 



and its boundary condition is 



(V$4) • (VGo) - ^(n X i)fig ■ (VGs) 



= 0. 



The equation for determination of the figure is 

Ae4 = -716^-^64, 
and its boundary conditions at the stellar surface are 



C4 



+ ^P3(sin^cos^) = -|^54 + E p^^r, 

^ l,m ''I 



5e 



= 0) = 0, 



The regularity condition at the center of the star is 

804 

and the normalization of G4 at the center of the star is 

e4(^ = 0) = 0, 

because we take ©o(^ = 0) = 1. 



6. 5th order 



The equation for the velocity potential is 



Go 



(V4>5) • (VGo) - -{n X • (VG4) 

?1 



and its boundary condition is 



(V$5) • (VGo) - ^(n X ^)fig ■ (VG4) 



= 0. 



The equation for determination of the figure is 

AGs = -nG^^Gs, 
and its boundary conditions at the stellar surface are 

- TT^^4(sin^cos^) = -^S, + '^Yr, 

^1 l,m ?1 

d^Go , , 5G5 , n ^ X 2ko , 



■(Ci 



^^P4(sin^cos^) = 'fs, - + l)|tll 



(^ = 0) = 0, 



The regularity condition at the center of the star is 

aGs 

and the normalization of G5 at the center of the star is 

G5(^ = 0) = 0, 

because we take Go(^ = 0) = 1. 
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7. 6th order 



n 



The equation for the velocity potential is 

A|>6 

and its boundary condition is 

(V$6) • (VGo) 

The equation for determination of the figure is 



(v$6) • (VGo) - -{n X ^)fig ■ (ves) 



-{n X Ofig ■ (ves) 



and its boundary conditions at the stellar surface are 



C6 



fi6 



2/x 



-'S'3-P2(sin6'cos(p) 



2Mtotal 2(l+p)2 ^1(1 +p) 

= 72i7--'^6j -2^(« + i)-i+^53y; +2^-1+1 

^ 



6 S<P 



(A58) 



(A59) 



(A60) 



1+p 

ym 



6 5v 



^3^2 (sin 9 COS if) + — r-Ps (sin 1 



'J q2 ^ ^ 

f ife(ei) = -^f ^ - 2^6) + + + 2)5S'53yr - E(' + 



1=0 



z=o 



a(l+p)^ 

Ha 



where 



5 = 



9 [fii ^ f^\^±n 



Ml Voq/ M2 



The regularity condition at the center of the star is 



(? = 0) = 0, 



and the normalization of Ge at the center of the star is 

e6(^ = 0) = 0, 

because we take ©o(^ = 0) = 1. 



(A61) 

' cos y) 
(A62) 

(A63) 

(A64) 
(A65) 



APPENDIX B: SOLUTIONS AT ZEROTH ORDER 

Since the configuration of star 1 is spherical at 0th order, the equation for determination of the figure becomes 
Lane-Emden equation, 



1 d / 2<^Qo 



For example, we obtain analytic solutions. 



00 = 1 - y (for n = 0), 
60 = ^ (for n = 1), 



(Bl) 

(B2) 
(B3) 
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or numerical solutions by solving the ordinary differential equation (Bl). 

The velocity field at 0th order is the same as the orbital motion so that the velocity potential is determined by 

V$o = ^R{^ X ^)orb- (B4) 

Then, the solution should be 

$0 = r sin 6* sin (y9. (B5) 

1 + p 

From this equation, we find that it is convenient to normalize $o ^ 

$0 = ^sin6'sini^. (B6) 
In the following, we use the normalized velocity potential defined by 



Here we show the useful relations: 



^ (B8) 



(V$o) • X = ,^ , (B9) 



V 

The orbital angular velocity can be calculated fi-om Eq. as 

->2 



1 

(1 ^V? 

( V$o) • (f^ X ^ ) = sin e cos (BIO) 



0./ + (BID 

J star 1 ^ R 



Therefore, we have 

i?3 



nl = ^, (B12) 



where Mtot = Mi + A'h- 

For 0th order, we have the internal and external gravitational potentials and their derivatives as 

C/z„t(0 = 0o + co, (B13) 

= ^, (B14) 

UeAO = J, (B15) 

= (B16) 



They become at the stellar surface as 



U„,t{^i)=c^, (B17) 
^(^i) = ^(6), (B18) 
UeA^i) = (B19) 



^(Ci) = -|. (B20) 
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Therefore, we can decide two constants at 0th order from the boundary conditions at the stellar surface as 

«o = -?f^(a), (B21) 

co = -a^(a). (B22) 



APPENDIX C: SOLUTIONS AT FIRST ORDER 
The equation for the velocity potential is 



A^i = -n 



1 



1 



(V$i) - X (CI) 



= -n(V$i) • (C2) 

where we have used the fact that V6o has only ^ component and {CI x ^) fig has only (p component. The boundary 
condition is 

(v$i)-(veo)|j^=o. (C3) 

When we expand $i as 

^i = Y,^'^ct>immr{e,^), (C4) 

l,m 

the boundary conditions are that ^^Vim is regular at ^ = and 



(6) = o. (C5) 



We cannot determine the absolute value of $i from the equation and the boundary condition. Therefore, $i should 
be zero. Note that $i = constant is a solution. However, since the physical value is a velocity field V$i, the constant 
solution is meaningless. 

Next, we consider the deformation of the figure at 0{e). The equation is written as 

AGi = -nej^-^ei. (C6) 

We expand 0i as 

ei = ^«V/™(0>l"(^,<p)- (C7) 

l,m 

The surface equation up to 0(e) is written as 

E = ^i + eSi{e,<p), (C8) 

where 

The gravitational potentials and their derivatives up to 0(e) are 
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Uint 

Uext 
dUext 



^ fi{3 + 2p) 
- e^4^ 'Jr + Co + eci, 



2(l+p)2 



dQ 



'^0 I \^ 



l,m 



■ —-e^{l + l)^Yi 



They become at the stellar surface as 



Uint{^) = Co + e 



ci 



m(3 + 2p) 



5i 



2(l+p)2 



l,m ^1 



5^ 



Then, the boundary conditions for 0(e) become 

/i(3 + 2p) _ _Ko 
2(1+P)2" ^? 



Cl 



5^ 



and also 6i should be satisfied with the conditions at the center of the star, 

t(o)^o, 
ei(o) = 0. 



(CIO) 

(Cll) 
(C12) 

(C13) 



(C14) 

(C15) 
(C16) 

(C17) 

(C18) 
(C19) 



(C20) 
(C21) 



1. The case of n = 

In this case, since we can express '•^Vfm = '"^^^im^'' by using a constant '•^•'Aj^, the boundary conditions can be 
written as 



Therefore, we obtain 



Cl - 



m(3 + 2p) 



2(l+p)2 

?1 



_ /x(3 + 2p) 
2(l+p)2' 

K.l:lm = 0, 

Gi = 0. 



(C22) 
(C23) 

(C24) 

(C25) 
(C26) 
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2. The case of n > 



In this case, we have the relation; 



(a) = -ir^i^i), (C27) 



from the Lane-Emden equation at the stellar surface. Then, the boundary conditions can be written as 

m(3 + 2p) ^ 



l.m, ^1 



Therefore, we obtain 



^^^^ ^ ■ 'e.' 



_ /x(3 + 2p) 

- 2(l+p)2' ^^^0) 

/^i:;™ = 0, (C31) 

Oi = 0. (C32) 



APPENDIX D: SOLUTIONS AT SECOND ORDER 

The equation for the velocity potential is 

A$2 = -n(Vl>2) • (Dl) 
Bo 

The boundary condition is 

(V$2)-(Veo)|j^=0. (D2) 

Then, $2 should be zero as in the case of $1. 

Next, we consider the deformation of the figure at O(e^). The equation is written as 

AG2 = -n@l-^Q2. (D3) 

We expand 02 as 

The surface equation up to O(e^) is written as 

S = a+e'52(^,v), (D5) 

where 

c, Q2(Cl) fy.f.-. 

The gravitational potentials and their derivatives up to O(e^) are 
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Uint = O + Co + e^C2, 

dUint de 



Uext 
dUext 



^ l,m ^ 



They become at the stellar surface as 



Uinti^) = Co + e^C2, 
dUint dOo X , 2 



^1 ^1 /.rr, ?1 



Then, the boundary conditions for 0{e'^) become 



C2 



-2i2 + 2^-^r; , 

^1 (,m ?1 



d^Go . . 962 , . _ 2ko Y^^, , T^ /t2:;m ^. 



l,m ^1 

and also 62 should be satisfied with the conditions at the center of the star, 

t(o, = o, 

62(0) =0. 



(D7) 
(D8) 

(D9) 

(DIO) 

(Dll) 
(D12) 

(D13) 
(D14) 

(D15) 
(D16) 



(D17) 

(D18) 



1. The case of n = 



In this case, since we can express ^^Vim = ^^^^imC' by using a constant ^"^^Aim, the boundary conditions can be 
written as 



C2 = j:{-^Mme. + ^)Yr, 



= Y^\{3-l)(^)Airn^[-^-{l + l) 



tl+2 
?1 



(D19) 
(D20) 



Therefore, we obtain 



C2 = 0, 
1^2:17X1 = 0, 

62 = 0. 



(D21) 
(D22) 
(D23) 
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2. The case of n > 



In this case, the boundary conditions can be written as 



Therefore, we obtain 



Ira 



(6)-(^ + i) 



1+2 



C2 = 0, 
Ii2:lm = 0, 

62 = 0. 



(D24) 

(D25) 



(D26) 
(D27) 
(D28) 



APPENDIX E: SOLUTIONS AT THIRD ORDER 



The equation for the velocity potential is 



A$3 = -n(V$3) 



en ■ 



The boundary condition is 

(v$3)-(veo)|j^=o. 

Then, $3 should be zero as in the cases of l>i and ^2- 

Next, we consider the deformation of the figure at O(e^). The equation is written as 



We expand ©3 as 



AOs = -ne^^Oa. 



e3 = Y^^'^immr{o,^)- 



The surface equation up to O(e^) is written as 



where 



S = a+e'S3(e,^), 
^ 93(6) 

' eo,j(a)' 

The gravitational potentials and their derivatives up to 0{e^) are 

Uint = Q + Co + 



C3 - 



2^ 



^ ^2p2(sin6'cos(^) 



^P2 (sin 6* cos ip), 



ext 



Yr 



l.m 



dUext _ Ko ,3V^n I T\ '^3:lm, ^rn 
l,m ^ 



(El) 
(E2) 

(E3) 

(E4) 

(E5) 
(E6) 

(E7) 

(E8) 
(E9) 

(ElO) 
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They become at the stellar surface as 



Uint{^) = Co + e' 



C3 - 



l+V 



P2 (sin 6 cos <p) 



2/x 



P2 (sin ^ cos if) 



Then, the boundary conditions for O(e^) become 



C3--—P2(sm0cos^) = --2-53 +2^-^11 . 



-(6) 



5$ 



■(cl) 



^^P.(sin^cos^) = -Ea + DlS^^^-, 



and also 63 should be satisfied with the conditions at the center of the star, 

63(0) = 0. 



(Ell) 
(E12) 
(E13) 

(E14) 

(E15) 
(E16) 



(E17) 
(E18) 



1. The case of n = 

In this case, since we can express '•'^Vim = ^^^^(mC' by using a constant the boundary conditions can be 

written as 



^3 - ^P2(sin0cos^) = + tS^) V, 

-77rT;A^2(sinecos^) = - ^r^M^t,-^ - (/ + 1)5 



;+2 



y-r, 



Therefore, we obtain 



where 



C3 = 0, 

63 = (3U2^^P2 (sin 6» cosy.). 



5/i 



2(i+p)er 



2(1+P)^ 



(E19) 
(E20) 



(E21) 
(E22) 



(E23) 
(E24) 



2. The case of n = 1 



In this case, since we can express '^V/m = ^^^Bimji by using a constant ^^^Bim and a spherical Bessel function ji, 
the boundary conditions can be written as 



26 
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APPENDIX F: SOLUTIONS AT FOURTH ORDER 



The equation for the velocity potential is 

n 



(v$4) • (veo) --{nx ^)fig ■ (vea) 



eo 



(V$4) • (VGo) + —(3^2^1(008^) Sin 2^' 



Therefore, we have the form of $4 and the equation for it as 

1.4 = (%(^)P|(cos6i)sin2(^, 



1 d 



e 



6 n dOo d 



2600 



(3^2. 



d^V^ d£,) $2 ^ ^ ^ 
The boundary condition is 

Next, we consider the deformation of the figure at 0{e^). The equation is written as 

Ae4 = -nQ^-^Qi. 

We expand ©4 as 



e4 = ^(4V;m(05^'"(^,v). 



The surface equation up to O(e^) is written as 



where 



„ _ 64(6) 



The gravitational potentials and their derivatives up to 0{e'^) are 



Uint 

Uext 
dUext 



P2 (sin 9 cos ip) + e 



C4 



(99 



.4 ,-2 



^^Ps (sin 6* cos ip) 



^^P3(sin^cos if), 



Ko 3 3K3:2n^_._rt _^ 4 



P2(sin0cos^)-e4^(; + l)5g^y, 



They become at the stellar surface as 



Uinti^) = Co - e^— ^^P2(sin6icos<^) + 
1+ p 



C4 + 



^1+2 

— P3 (sin 6 cos 



dUjnt 



dBo 



(6) + e= 



Sa—^i^i) + -Er^(6) - ^ , , P2(sm6>cos(p) 



^4 



de 

S: 



(6) 



de 
964 



(6) 



9^ ' 6(1 



jP3(sin^cos 

£/e.*(-) = f{l-e'^f- e^|i) + e^^P2(sin^cos^) + e^E Tm^^ ' 



(") = -72(l-e^ -e ^) -e''^4--P2(sin0cos^)-e 2^{l + -f+^^i 
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Then, the boundary conditions for O(e^) become 



C4 



+ ^P3(sin^cos^) = + E ^yr^ 



Ps{.meco.^) = ^-^S,-J2{l + l)^Yr 

^1 l.rr,. ^1 



6(1 +P) 



and also 64 should be satisfied with the conditions at the center of the star, 

64(0) = 0. 



(F18) 
(F19) 



(F20) 
(F21) 



1. The case of n = 

In this case, since we can express ^^Vim = ^^^^im^'' by using a constant ^^^Aim, the boundary conditions can be 
written as 



^4 + Y^P3(sin^cos^) = + ^)Yr, 

^ l,m ^1 

-^P3(sin^cos^) = - - {I + 1)5 



Im 



yr, 



Therefore, we obtain 



where 



C4 = 0, 

64 = ^^^Aj,i^Po,{^uiQzo^^>), 



K4:3 = - 



7m 



4(i+p)er 

4(1 



(F22) 
(F23) 



(F24) 
(F25) 



(F26) 
(F27) 



2. The case of n = 1 

In this case, since we can express '^Vim = ^'^^Bimji by using a constant ^^^Bim, the boundary conditions can be 
written as 



C4 + -^P3(sin^C0S^) = J^\-^Mmjl{^l) + ^ 

l,m 

^-^P3(sin^cos^) = E[-^'^^''"|(^i) - + 



, '^4:in 



Therefore, we obtain 



C4 = 0, 

e4 = (^)B3j3(0^3(sin^cos^), 



(F28) 
(F29) 



(F30) 
(F31) 



where 
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K4:3 = - 



3(1 



3(1 +p) 



(21 - ), 



jsiO - ^ (15 - 6^') sin^ - ^(15 - e) cosC 



3. The case of n > 



In this case, the boundary conditions can be written as 

C4 + P3 (sine cos ^) = V[-(4Vj„(a) + 5 



hn 



Y, 



3/1 



6(1 +P) 



P3(sin6'cos(/5) = 



Im 



(ei)-(^ + i) 



I 1 



1+2 



Therefore, we obtain 

C4 = 0, 

64 = W^3(0^3(sinecos<^) 
and the ordinary differential equation for ^^^i/js is written as 

1 d d\ 12 



The coefficient of ^^Vs and the constant k;4:3 are determined from two equations; 



d^ 
d^ 



(6) = - 



K4:3 = 



('V3(6) + Y 



4. Octupole moments 

At this order, we can obtain the octupole moments, but we do not represent them here. 



APPENDIX G: SOLUTIONS AT FIFTH ORDER 



The equation for the velocity potential is 

A$5 = 



n 

n 



(v$5) ■ (voo) - ^(n X i)fig ■ (ve4) 

(Vl>5) • (VGo) - ^(4)^3(^3 (cos0) sinvJ- ip|(cos0) sin3^j 



Therefore, we have the form of $5 and the equation for it as 

|>5 = (5)^3 (^) P3i(cos6l)sin</? - ip3^(cos6l)sin3(^ 



1 d / 2 d \ 12 n dQo d 
'^d^y^ 'd£,)~^^ Q^~dfdi 



(5^3 



4^60 



(4)^3. 
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The boundary condition is 



Next, we consider the deformation of the figure at 0{e^). The equation is written as 



We expand 65 as 



AOs = -ne^-^Gs 



The surface equation up to 0{e^) is written as 



where 



S = 6 + e^Saie, ^) + e^S4{e, if) + e^S5{e, if), 

^ ©5(6) 

or, = — - 



The gravitational potentials and their derivatives up to O(e^) are 



^2p2(sin6lcos(/7) + 



^^P3(sin^cos ip) 



99 



C5 



^'*P4(sin0cos (fi) 



—e 



e^^TT^ — rC-P2(sin^cos(^) + e 



4 '^M i-2 



^^P3 (sin 6 cos ip) 



^ ^3p4(sin6'cos(p), 



Uext 
dUext 



eAi+p) 

1^0 . s '''3:2 



^ + e -^^^(sin 61 cos v>) + e*-^P3(sin 6» cos v-) + e'^ 2^ -^j^y; 



KO 3 3K3:; 



P2(sin^cos^) - e^^^P3(sin0cos^) - ^(Z + l)^^!^'". 

^ l,m ^ 



They become at the stellar surface as 

Uint (S) = Co - P2 (sin 6 cos + — ^ P3 (sin 6» cos + e' 

1 + P I + p 



dQo 
d^ 



(6) + ' 



C5 



1 +p 



P4(sin0cos 



2/^ 



'"^^ 6(1 +P) 

3/i 



P2(sin^cos lyj) 



P3(sin^cos 



^.^(a) + ^(6)-^^P4(sin.cos^) 



3 'S'3 



S'4 



^% h-O /, 3»33 4J4 5'^5\ , 3«3:2r) / . ^ n , 4 K4:3 „ / . ^ x , 5 V"^ 

Uexti^) = ~ ~ TtJ +e^^2(sm6'cos(/p) + e*^P3(sm6'cosv?) + e''2^ 



dU, 



ext /I— 1 



(H, = -|( 



6 6 a 



2^ g3 2*^3 ^4 25*4 ^5 



6 



a 



;3k3:; 



■P2 (sin ^ cost/?) — -^|^P3 (sin ^ cosy) 



Lm. ?1 
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Then, the boundary conditions for O(e^) become 

— /ill Sin 6/ COS = — -0-D5 + 
r> fr 



C5 - — -P4(sin^cos^) = --5-55 + 2^ -i+r^; , 



and also 65 should be satisfied with the conditions at the center of the star, 

65(0) = 0. 



(G18) 
(G19) 



(G20) 
(G21) 



1. The case of n = 

In this case, since we can express ^^V/m = ^^^^;mC' by using a constant ^^^Aim, the boundary conditions can be 
written as 



^5 - YT1^^4(sin^cos^) = ^(-(^U^^^i + p^)Yr, 



l,m 



Y7 



Therefore, we obtain 



where 



C5-O, 

05 = (5U4^''P4(sin6»cos<^), 



(5)^4 = 



2(i+p)er 

2(1 



(G22) 
(G23) 



(G24) 
(G25) 



(G26) 
(G27) 



2. The case of n = 1 

In this case, since we can express '^Vim = '"^^Bimji by using a constant ^^^Bi^, the boundary conditions can be 
written as 

(G28) 
(G29) 



Therefore, we obtain 



C5 = 0, 

65 = (5)B4j4(0A(sin^cos^), 



(G30) 
(G31) 



where 
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K5:4 



(l+p)(15-a' 



(15 -e? 



945 \ 



i4(0 = [(105 - 45^' + sinC - e(105 - 10^^) cos^ 



3. The case of n > 



In this case, the boundary conditions can be written as 

C5 - :r^P4(sin^C0S^) = ^[-(^^(^(a) + ^ 



Im 



Y, 



4/i 



ei(i+p) 



P4(sin0cos^) = ^[-^l^fe) - + 1)1^^^^ 



Yr 



Therefore, we obtain 



C5 = 0, 

05 = (5)^4 (0^4 (sin 6* cos <^), 
and the ordinary differential equation for ^^^'^/^ is written as 

1 d 



The coefficient of '•^■'V'4 the constant K'^-a arc determined from two equations; 



5(^V4(a)+6^(a^- 



«;5:4 = Cl 



4m 
1+p 



(a) 



4. Orbital angular velocity 



At this order, the orbital angular velocity has the 5th higher order term than the monopole term, 
velocity is 



0^ = 



i?3 



1 



2B? \ Ml 



Mn 



For the calculational convenience, we rewrite the above equation as 

Mtot, 



ft' 



i?3 



where 



S = 



2al 



Ml Vao/ M2 



(G32) 
(G33) 
(G34) 

(G35) 
(G36) 

(G37) 
(G38) 

(G39) 

(G40) 
(G41) 

The total angular 
(G42) 

(G43) 

(G44) 
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APPENDIX H: SOLUTIONS AT SIXTH ORDER 

The equation for the velocity potential is 



©0 

n 



(v$6) • (VGo) - ^{n X • (ves) 



(V$6) ■ (VOo) - Y^(''^V4{p|(cose) sin 2^ - ^Pi{cos9) sin4(^}] . 
Therefore, we have the form of $6 and the equation for it as 



|>6 = (6)^4(0 [^4 (cos 6») sin 2ip - -P4 (cos 9) sin 4v? 
1 d d\ 20 n dGo d 



i d / 2 d \ 

£2 W~e 



ed^y' d^J e ©0 rf£ 



12a©o 



(5V4. 



The boundary condition is 



Next, we consider the deformation of the figure at O(e^). The equation is written as 



We expand ©e as 



A06 = -n©"-i©6 - -n{n - 1)©^^©^ + ^A(%(0f|(cos6') cos2(^. 
2 6 



The surface equation up to 0(e®) is written as 

S = 6 + e^Ssie, ^) + €^54(0, ^) + e'S^ie, ^) + £656(0, 

where 

1 



56 = - 



©o,c(6) 



^^(£1) + 53^(6) + ©6(a) 



2 



The gravitational potentials and their derivatives up to 0(e®) are 



Uint = © + Co - e 



3 ^ f P2 (sin cosy) + 6 ^3 ^, g^P3(sin^cosy) 



e?(l+f) 



e?(l+p) 



^'|4(r^^'^4(sin0cosy) 



C6 - 



3mJ^ii /do 



(-)'- 

van/ 



/x5 



2M..agV.oy 2(1+,)^ .^£^P.(sin.cos.)-|(%P|(cos.)cos2, 



3 ^ £P2(sin0cos^) + e ^3,^' ^. g^PaCsingcos^) - 



e?(i+P) 



ef(l+p) 



£'^P4(sin0cos(p) 



e*P5(sin e cos <p) - ^ ^^P|(cos ^) cos 2^^' 



Uext 
dUext 



_ ^ e^:^P2(sin0cos^) + e*:^P3(sin0cos^) + e^:|^P4(sin0cosv5) + 5] ^rr, 



3 3«;3:2„.. ^ X 4^'^4::3 jj , . n \ 5 ^'''5:4 „ / . n \ 

--^ — e P2 (sm 6^ cos y ) — e P3 (sin cos y?) — e P4 (sin cos y) 



-e6B^ + l)5gil^^ 



(HI) 
(H2) 

(H3) 

(H4) 

(H5) 

(H6) 
(H7) 

(H8) 
(H9) 



(HIO) 

(Hll) 

(H12) 

(H13) 
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They become at the stellar surface as 



UintC^) = Co - e 



3 M 



l + p 



2S: 



1 + ^e^]P2{sin9cos<f) + e 



4 M 



l+P 



' \ 3 



C6 



Van/ 



/i(5 



P3 (sin cos (fi) — e 



5 



P4(sin0cos(p) 



on ^ ^2 + 7^^^5(8111 e cos - ^(%(a)P|(cose)cos2vp 



2Mtota^ Vao 

-^"^ = + ^ [^^^(^^) + - a(TT^(^ + ^e3)p.(s.n^cos^) 

3^ 



(H14) 



54^(6) + ^(a) 



9^ 



6(1 +P) 
4/x 



90 



P3(sin^cos 
P4(sin^cos<p) 



(6) + 



lc,2rf'©o 



(ei) + 53^(6) + 



dee 



(a) 



5/i 



2m ^(4)02 



(6)P|(cos6i) cos2<^ 



a ^ a J J 



a(i+p) 



5 353\ K3: 



P5 (sin d cos 



1 - ^ ^P2 (sin 9 cos (^) 



a ^ a^ 



+e^^P3(sinecos^) + e5^P4(sin^cos^) + ^ 



5^5:4 



1 — f " 

,1 ^ a a 

44«;4:3 D / . /) X 

-e — r^P3(smt'cos 1^) — e 



^1 l,m ^1 

IHaJ -'irlJ I'"' — J^P2(sm0cos^) 



(H15) 



(H16) 



(H17) 



Then, the boundary conditions for O(e^) become 

C6 



C (i^lp) -^3^2 (sing cosy)) + j^P5(sin 6* cos (p) - |^(%(6)P|(cos 6') cos2(^ 



= 72 ( ;r - '^6 j - — 3-53P2(smecos<p) + 2^ -j^l^ 



5293 



2n 

"a 



(H18) 

P5 (sin 6 cos 

(ei)P^(cos(^) cos2^ = -|(^ - 2^6) + i^53P2(sin0cos^) - Y.{1 + l)^Fr, (H19) 



^2 -(a) + 2^3^(a) + ^3^(a) + ^(a) - 



53P2 (sing cos ^)+ 



a(i+p)' 



a^ 



where 



C6 = C6 



2Mtotal\ao) 2{1 +pf' 
Also Oe should be satisfied with the conditions at the center of the star, 

06(0) = 0. 

In the calculation of this section, we need the expression of [P2(singcos(^)]^. It is written as 

2 18 2 1 

rP2 (sin 6 cos (p)] = — P4 (sin 9 cos </?) + - P2 (sin 9 cos t^) + - . 

35 7 5 



(H20) 

(H21) 
(H22) 

(H23) 



35 



1. The case of n = 



In this case, since we can express ^^^2 = ^^^^2^^ and ^^^im = ^^^^im^'i by using the constants '^^^2 and ^^^Aim, 
the boundary conditions can be written as 

C6 + ---^P5(sin6'cos(/p) - ^(%(a)P|(cos6l) cos2v? 

+ 5;(-CU,„{! + ^5±i)yr, (H24) 

— (sin e cos ^) - 1^ ^ (6 )P| (cos ^) cos 2^ 

= rWTpj ^ 2''^' - (sS^^^^^^^'^^^^) + -P2(sm.cos^) + -) 

+ ^ [(3 - -{1 + ] Yr. (H25) 

l,m ^1 

Therefore, we obtain 

Ge = ('^U2^2P2(sin6'cos(p) + ('3U22^^P|(cos e*) cos2ip + ^^U^^'^Piisine costp) 

(sin ^ cosy), (H27) 



where 



'"^5 = (H31) 



8(1+ rt£ 



2. The case of n = 1 



In this case, since we can express ^^^2 = ^^^P2j2 and ^^^ipim = ^^^Bi^ except for ^^^22 by using the constants ^^^52 
and ^^^Bim, the boundary conditions can be written as 



C6 + — ^Pg (sin 6*008 i^) - ^ (4)02 (^^)p| (cos 6*) 008 2^5 
1 + P Ci 



0o,«(a) 



2/i 



dj2 ^ , 3k3:2 



6(1 +P) it 



(18 2 1 \ 

— P4 (sin ^ cost/?) + -P2 (sin cost/?) + - ) 



6(1 +P) 



P5(sin0cos <p) 



2/x 



(6)P2^(cos6i) cos2(/7 



0o,«(6) L e?(l+p) 



+ 



(H32) 
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/ 18 2 1 \ 

X ( —P4 (sin ^ cose/?) + - P2 (sin ^ cose/?) + 7) 



SI 



Therefore, we obtain 

C6 = 



96 = ^^'^B2j2{0P2 (sin e cos ip) + ^^^22(0^^ (cos 6») COS 2(p + (6)B4i4(0'P4(sin 6» cos cp) 
+^®^S5j5(e)^'5(sin0cos¥'), 



wliere 



405/^2 



7(1+^)2(15-^2)' 



(l+p)(105-10Ci2) 



2\' 



j5(^) = — (945 - 420^2 _^ 15^4) gj^^ _ ^(-945 _ -^Q5^2 _^ ^^^^ 



and the coefficient of '^^22 is determined from the equation; 



3(^2(6)4-6^(6) = ^ 
"? 6 



3(%(6)+6^(6) 



3. The case of n > 



In this case, the boundary conditions can be written as 



C6 



2/x. 



-P5(sin6'cos(p) - -f-(%(6)-P|(cos6l) cos2v? 
1 + P 6 



('¥2(6) 



2/t 



00, ?(6 

1. m SI 



6(1 +P) 



(6) + (^i^P4(sin0cos v) + |P2(sin^cos</7) + 



Y, 



I > 



P5 (sin e cos (p) - ^ (6 )P2 (cos 0) cos 2(^ 



6(1 +P) 



6 rfC 

2^ / 3 



1 



(^^2(6) 

0o.c(6) L eAi+p) ' V6' 2eo,c(6) 



(6))('V2(6) 



+ 



(i^2 



(6)- 



12k 



3:2 



6' 



+ 



X ('i^P4(sin0cos(p) + ^P2(sin0cos(p) + 
\35 7 5 



l.m 

Therefore, we obtain 



^(6)-(« + i) 



?1 



Yr. 



6 



"0o"'(6)(('V2(6)) - ('Vo(6) - 6-ir^(6), 



1060,4(6) """'^ rfC 

96 = (%(0 + (0-^2 (sin 6* cos (^) + (6V22(O^I(cos6')cos2(^+ (6)^4(0-P4(sin6'cosv5) 
+(6V5(0'P5(sin^cos^), 



37 



where the forms of ^^Vi ^-I'S determined by the equations 



i a / 2 a \ 






n(n 

10 ^ 


-1)9 


^'((^¥2)^ 


r ^ rl / rl 
( ^2 

-ed^y d4 






1(6)V.2 - 




-i)er2pV2f, 


r T rl / rl 
I r2 

-ed^y d4 








9// 


r 1 W / /7 \ 


" 1 d z 2 d 
-ed^y d4 


\ 20- 


(6)^74 = -ne^ 




35 ^ 


r^-l)e^^((3V2)^ 


" 1 d / 2 d 
-^d^y d^ 


\ 30- 











and the coefficients of ^^"^ipi are determined by the boundary conditions; 



7eo,aei) 

2/i 



3(«V22(6)+6^(6)-, 



3(%(a)+a^(6) 



3590,^6) 
11^ 



(H45) 
(H46) 
(H47) 
(H48) 
(H49) 

(H50) 
(H51) 
(H52) 
(H53) 



APPENDIX I: THE CHANGE IN THE CENTRAL DENSITY 

In this section, we present the expression of the change in the central density at O(e^) by using the first tensor 
virial relation. 



First, we divide the first tensor virial relation ( 2.28 ) 



n n 



^Hi + ^n2 + (PFsei/)tot + {W^^t)tot + mot = 0, 



(11) 



into three parts of 0(e ), 0(e) and 0(e ). After that we demand that these three parts should be zero independently, 
(a) The 0th order of e: 
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(12) 



If we assume that both of the terms concerned with star 1 and star 2 become zero independently, we have two relations; 

5-n 1 rii 



l + «e?(0o,c(a))2 Jo 
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2ol+n 



1, 



1 + "' ei'(eo,«'(ei))'^o 

The above relations are proved using the Lane-Emden equations 

(b) The 1st order of e: 

It is clear that this term becomes zero. 

(c) The 6th order of e: 
By using the relations; 



(13) 
(14) 
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C6 = 



Pc 

3 -n'\ ^, - ,^,.Sp'^ 



2n 
l+p 
P 



eieo,e(6), 



A* =-(i+p)ei0o,e(e;), 

we can express the equation at order as 
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If we assume that both of the terms concerned with star 1 and star 2 become zero independently, we obtain two 
relations; 
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where we have used Eqs. (ttS^ and 
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Tables 



TABLE I. The functions and their derivatives to represent the equiUbrium figure for n = 0.5. (I) 



1) = ()..) 


c 
s 










ip.i 


a ipz/ ^Z, 




Uj {^4 / 


n nn 
u.uu 


i.UUU 


u.uuu 


n nnn 
u.uuu 


n nnn 
u.uuu 


000 
u.uuu 


000 

u.uuu 


000 

u.uuu 


000 
U.UUU 


10 

U. J-U 






o.o 1 i-y OJ 


n 1 074 

U. ±U ( ^ 


1 HI fir 4"! 


Q q4fi(' "Xx 


4:.UUO\^ U j 


1 fini ( d^ 

± .UU J- \^ 4: f 






-fi fi^l'^r p^l 


9 1 Ap./ 9"! 


0.2143 


-1 n^9('-H"i 


-1 'i7fil' 9"! 


fi 4no(' K\ 

U.4:UU\^ j 




^0 






4 89n(' 9") 


n ^90^ 

u.ozuo 


o.o^^y J 


Q Krinf 9~\ 


O.ZOU\^-4: ) 


4 '^1 1 ^ "i"! 


n 40 


n Q7'id 
u.y 1 04: 


-u. ±ozo 


a ^AS,( 9"! 


n 494Q 


-o.oou l^-oj 


R 970^' 9^ 


1 091 (' '^'l 

1 .UZ ± \^-OJ 


1 01 9^ 
± .uiy i^-z^ 


^0 
u.ou 


u.yoou 


-U. lU^U 


n 1 ^^1 


n ^977 


1 fi'^4(' 9"! 


Q 7'^'^^' 9~1 

y. 1 ooi z J 


9 488 f '^'l 

z.^oo^^-o ) 


1 Q8'^l' 9^ 
± . you\^-z J 


n Ro 

u.uu 




n 1 Q^^4 


n 1 QHQ 

u.iyuy 


U.OZoU 


9 81 4/' 9^ 


n 1 ^Q7 


^ 1 4fi^' 




n 70 

U. 1 u 


U.£7l£70 


n 997fi 


n 9^18^ 


n 79^^ 


4 4^9^ 9"\ 

— ft.^OZl -Z J 


n i8Qn 

-u. J-Oyu 


y.ouu ) 


U.^UUl^-Z J 


80 
u.ou 


U.OcJOU 


n 9^81 


u.oooy 


n 81 Q7 
u.oiy 1 


fi fil 7f 9"! 


n 94 


1 filfil' 9") 


8 01 7(' 9~> 
o.ui ( j 


n Qo 
u.yu 


u.ou ( o 


-U.Zo / / 


n /1 99/1 


u.yuyy 


Q '^7A ( 9^ 

-y .0 / ^i^-z J 


"^07^ 
-U.OU / 


9 t=;77f' 9"^ 


1 1 "^4 


1 00 

1 -UU 


u.oo / o 


U.OlUtJ 


n 111 78 

U.OJ. 1 


u.yyou 


n 1 97Q 
-u. i-Z / y 


'^7^0 
-U.o ( oy 


O.yj-Ui z 1 


1 ^4^ 
u. 10^0 


1 10 

J.. ±u 




-0.3441 


0.6214 


1 n7fi 


_n 1 fiQI 


-fl 44QS 




n 9040 


1 90 
1 .zu 


U. / uo / 


-U.o / u / 


n 7^98 


1 1 c^l 
X. 101 


n 91 80 
-U.ZloU 


p;987 


8 090f' 9^1 
o.UZUl^-Z j 


9fi94 
u.zoz^ 


1 '^0 
1 .ou 






n 8'ii 4 


1 91 Q 
± .z±y 


n 97'i0 
-u.z J ou 


fil 1 8 
-u.ui 10 


n 1 nQ7 

u. luy ( 


u.oouz 


1 40 
± .^u 


U.DoyO 


0/11 Q7 


n 07fi/i 
u.y / 


1 980 
i.ZoU 


-U.O^UO 


fiQ8^ 

-u.uyoo 


n 1 /Ififi 
u. 1^00 


/1 078 
U.^U / 


1 "^0 

J. .UU 


n fl4.fi 


-0.4420 


1 107 


± .000 


-0.4148 


_n 787Q 
u. i ( y 


1 QI fi 
u. J.y j-u 




1 fiO 

1 .UU 






1 94^ 


1 ^77 
1.0/ / 


n 4Q81 


n 878Q 
-u.o 1 oy 


n 94^Q 
u.z^oy 


u.oyzu 


i. / u 


U.004U 


n /I Q 1 /I 
-U.4oi4 


1 QC9 


1 /111 


n Knot; 
-U.oyUD 


-U.y ( Uo 


U.oiUo 


U.OyyO 


1.80 


0.5051 


-0.4982 


1.525 


1.433 


-0.6922 


-1.061 


0.3862 


0.8158 


1.90 


0.4545 


-0.5128 


1.668 


1.442 


-0.8027 


-1.149 


0.4739 


0.9404 


2.00 


0.4026 


-0.5250 


1.813 


1.436 


-0.9219 


-1.233 


0.5745 


1.072 


2.10 


0.3496 


-0.5346 


1.955 


1.412 


-1.049 


-1.310 


0.6885 


1.209 


2.20 


0.2958 


-0.5413 


2.094 


1.367 


-1.183 


-1.376 


0.8164 


1.349 


2.30 


0.2414 


-0.5448 


2.228 


1.296 


-1.324 


-1.427 


0.9582 


1.487 


2.40 


0.1869 


-0.5447 


2.352 


1.190 


-1.468 


-1.456 


1.114 


1.617 


2.50 


0.1326 


-0.5403 


2.464 


1.037 


-1.614 


-1.451 


1.281 


1.728 


2.60 


7.904(-2) 


-0.5306 


2.557 


0.8061 


-1.756 


-1.389 


1.458 


1.798 


2.70 


2.674(-2) 


-0.5139 


2.620 


0.3959 


-1.887 


-1.189 


1.637 


1.753 


2.75270 


0.000 


-0.5000 


2.628 


-0.3645 


-1.942 


-0.6776 


1.724 


1.368 



41 



TABLE II. The functions and their derivatives to represent the equilibrium figure for n = 1. (I) 
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TABLE III. The functions and their derivatives to represent the equiUbrium figure for n = 1.5. (I) 
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yj . ij%jyjyj 


2.60 


n 2823 


-0.3313 


0.7229 


0.2654 


-0.3876 


-0 3365 

yj . Kjfjyj'j 


0.2463 


0.3241 


2.70 


0.2497 


-0.3204 


0.7488 


0.2519 


-0.4219 


-0.3485 


2802 


0.3531 


2.80 


0.2182 


-0 3088 


0.7733 


2393 


-0.4573 


-0 3606 


0.3170 


n 3836 




u.io / y 


-u.zyo 1 


u. / yo / 


n 9978 


n 4Q4n 


-u.o 1 oU 


u.oooy 


n 41 1^7 


3.00 


0.1589 


-0.2843 


0.8189 


0.2180 


-0.5319 


-0.3860 


0.4002 


0.4497 


3.10 


0.1311 


-0.2715 


0.8403 


0.2102 


-0.5712 


-0.3999 


0.4470 


0.4859 


3.20 


0.1046 


-0.2588 


0.8611 


0.2047 


-0.6119 


-0.4151 


0.4975 


0.5246 


3.30 


7.931(-2) 


-0.2460 


0.8814 


0.2021 


-0.6543 


-0.4323 


0.5520 


0.5666 


3.40 


5.534(-2) 


-0.2334 


0.9016 


0.2031 


-0.6985 


-0.4519 


0.6109 


0.6124 


3.50 


3.262(-2) 


-0.2212 


0.9221 


0.2086 


-0.7448 


-0.4752 


0.6746 


0.6633 


3.60 


1.109(-2) 


-0.2094 


0.9435 


0.2206 


-0.7937 


-0.5041 


0.7438 


0.7213 


3.65375 


0.000 


-0.2033 


0.9557 


0.2318 


-0.8213 


-0.5240 


0.7835 


0.7575 
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TABLE IV. The functions and their derivatives to represent the equihbrium figure for n = 2. (I) 



n = 2.0 


f 
s 


©n 






Li iy2 1 






f 4 


d'^^^ibA Idf 


0.00 


1.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


n QQS3 




1 195^-31 


9 oiniJA 


-1.631(-5) 


-4.890(-4) 


3.016(-7) 


1 906 f- 5^1 


0.20 


0.9934 


-6 614f-2) 


4.481(-3) 


4 456f-2) 


-1.301(-4) 


-1.945(-3) 


4.812(-6) 


9 607f-5) 


0.30 


0.9851 


-Q K22(-'l\ 


1 001 ('-9^ 


6 589('-2'l 


-4.366(-4) 


-4.337(-3) 


2 425 ('-51 


3.221(-4) 


0.40 


0.9738 


-0.1292 


1 7(\')(-')\ 


8 61 3('-9'; 


-1.027(-3) 


-7.612(-3) 


7.617(-5) 


7.563(-4) 


0.50 


0.9594 


-0.1586 


2 71q/-2"1 


0.1050 


-1 986/-3"l 


-1.170(-2) 


1.845(-4) 


1.460(-3) 


0.60 


0.9421 


-0.1863 


3 85fif-9"l 


0.1222 


-3 3Q9f-3"l 


-1 653f-9'l 


3.789(-4) 


9 4S6/-3'l 


0.70 


n 9999 


-0.2119 


5 1 5fi('-9'l 


0.1375 


-5 31 3('-3'l 


-9 iqqi'-9') 


6.941(-4) 


3 883('-3'l 


0.80 


n SQQ8 


-0 9359 




0.1508 


-7.807(-3) 


-9 799l'-9') 


1.169(-3) 


5 6S9('-3') 


0.90 


0.8752 


-0.2561 


8 1 65f-2'l 


0.1620 


-1 093f-2'l 


-3.443(-2) 


1.847(-3) 


7.934(-3) 


1.00 


0.8487 


-0.2745 


9 833('-9"l 


0.1712 


-1.470(-2) 


-4.120(-2) 


2.771(-3) 


1.064(-2) 


1.10 


0.8204 


-0.2904 


0.1158 


0.1782 


-1.917(-2) 


-4.821(-2) 


3 QQ1 


1 3S3('-9') 

X . OOO I Zi J 


1.20 


0.7907 


-0 3036 


0.1339 


0.1833 


-9 435f-9'l 


-5 53Si'-9'; 


5 553f-3'l 


1.750(-2) 


1.30 


0.7598 


-0.3142 


0.1524 


0.1864 


-3 095('-9'l 


-6 969f-9~) 


7.507(-3) 


9 166('-9') 


1.40 


0.7279 


-0.3224 


0.1711 


0.1879 


-3 fi88f-9"> 


-fi 9S6/-9'1 


9 909 ('-3'! 


9 639l'-9') 

Z.VJOZil Zi ( 


1.50 


n fiQ54 


-0.3281 


0.1899 


0.1878 


-4.422(-2) 


-7.706(-2) 


1.279(-2) 


3.147(-2) 


1.60 


0.6624 


-0.3316 


9086 


0.1863 


-5 99S/'-9"l 


-8.417(-2) 


1.621(-2) 


3.710(-2) 


1.70 


0.6291 


-0 3330 


0.2272 


0.1838 


-6 1 05('-9'l 


-9.117(-2) 


9 099('-9'l 


4 S93('-9') 

Ti.OZiOV Zi 1 


1.80 


n 5958 


-0 3395 


0.2454 


0.1803 


-7.051(-2) 


-Q 8n5l'-9~) 


2.487(-2) 


4.984(-2) 


1.90 


0.5627 


-0 3309 


9639 


0.1762 


-8 0fifif-9"> 


-0.1048 


3 091 f-9"> 


5 693^-9'! 

lJ-\jiJKJ I ~Zl 1 


2.00 


0.5298 


-0 3963 


2806 


0.1715 


-9.147(-2) 


-0.1115 


3 69Sf-9^ 


6 453(^-2') 


2.10 


0.4975 


-0.3211 


0.2975 


0.1666 


-0.1029 


-0.1180 


4.313(-2) 


7 969('-9'; 


2.20 


0.4656 


-0.3147 


0.3139 


0.1616 


-0.1151 


-0.1245 


5 0S9('-91 


8.124(-2) 


2.30 


0.4345 


-0.3073 


3998 


0.1567 


-0.1279 


-0.1310 


IS Q40l'-9'l 


9 n39l'-9~l 

i/.UOi/ 1 Zi 1 


2.40 


0.4042 


-0.2991 


0.3452 


0.1519 


-0.1413 


-0.1375 


fi «q9/.9"\ 


0.1001 


2.50 


0.3747 


-fl 9Q09 


3609 


0.1475 


-0.1554 


-0.1441 


7.944(-2) 


0.1104 


2.60 


0.3462 


-fl 9S0S 


0.3748 


0.1435 


-0.1701 


-0.1509 


9 109('-9'l 


0.1214 


2.70 


0.3186 


-0.2710 


3889 


0.1401 


-0.1855 


-0.1578 


0.1037 


0.1330 


2.80 


n 9990 


-0.2610 


0.4028 


0.1372 


-0.2017 


-0.1650 


0.1177 


0.1454 


2.90 


0.2664 


-0 9508 


0.4164 


0.1351 


-0.2185 


-0.1724 


0.1328 


0.1586 


3.00 


0.2418 


-0.2406 


0.4298 


0.1337 


-0 9369 


-0.1803 


0.1494 


0.1726 


3.10 


0.2183 


-0 9305 


0.4432 


0.1330 


-0.2546 


-0.1886 


0.1674 


0.1875 


3.20 


0.1957 


-0.2204 


0.4565 


0.1332 


-0.2739 


-0.1973 


0.1869 


n 9035 


3.30 


0.1742 


-0.2106 


0.4698 


0.1341 


-0.2941 


-0.2067 


0.2081 


9905 


3.40 


0.1536 


-0.2010 


0.4833 


0.1359 


-0.3153 


-0.2166 


0.2311 


0.2387 


3.50 


0.1340 


-0.1917 


0.4970 


0.1386 


-0.3374 


-0.2273 


9559 


n 9589 

\j . zooz 


'i fin 


U.llOo 


n 1 897 


U.Olll 


n 1 A99 


-U.oOU ( 


-u.zoo i 




n 97Q1 


3.70 


9.742(-2) 


-0.1740 


0.5255 


0.1466 


-0.3852 


-0.2510 


0.3118 


0.3015 


3.80 


8.043(-2) 


-0.1658 


0.5404 


0.1519 


-0.4110 


-0.2641 


0.3431 


0.3255 


3.90 


6.425(-2) 


-0.1579 


0.5559 


0.1581 


-0.4381 


-0.2783 


0.3769 


0.3513 


4.00 


4.884(-2) 


-0.1504 


0.5721 


0.1653 


-0.4666 


-0.2935 


0.4134 


0.3790 


4.10 


3.416(-2) 


-0.1433 


0.5890 


0.1733 


-0.4968 


-0.3098 


0.4528 


0.4087 


4.20 


2.016(-2) 


-0.1367 


0.6068 


0.1824 


-0.5287 


-0.3274 


0.4952 


0.4407 


4.30 


6.811(-3) 


-0.1304 


0.6255 


0.1924 


-0.5623 


-0.3463 


0.5410 


0.4751 


4.35287 


0.000 


-0.1272 


0.6358 


0.1980 


-0.5809 


-0.3569 


0.5666 


0.4944 
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TABLE V. The functions and their derivatives to represent the equiUbrium figure for n = 0.5. (II) 



n = 0.5 



f 

s 






















0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


1.720(-10) 


1.032(-8) 


5.400(-3) 


0.1080 


9.847(-3) 


0.1969 


2.061 (-6) 


8.245(-5) 


-1.312(-7) 


-6.558(-6) 


0.20 


1.104(-8) 


3.317(-7) 


2.158(-2) 


0.2155 


3.935(-2) 


0.3929 


3.299(-5) 


6.601 (-4) 


-4.195(-6) 


-1.049(-4) 


0.30 


1.264(-7) 


2.535(-6) 


4.846(-2) 


0.3220 


8.838(-2) 


0.5873 


1.672(-4) 


2.231(-3) 


-3.183(-5) 


-5.301(-4) 


0.40 


7.150(-7) 


1.078(-5) 


8 594f-2) 


0.4272 


0.1567 


0.7792 


5 291f-4'l 


5 299f-3) 


-1 339f-4~) 


-1 672f-3'l 


0.50 


2.752(-6) 


3.330(-5) 


0.1338 


0.5305 


0.2441 


0.9678 


1.294(-3) 


1.038(-2) 


-4.080(-4) 


-4.072(-3) 


0.60 


8.309(-6) 


8.410(-5) 


0.1920 


0.6315 


0.3502 


1.152 


2.689(-3) 


1.799(-2) 


-1.013(-3) 


-8.419(-3) 


0.70 


2 123('-5'l 


1.850(-4) 


0.2601 


0.7297 


0.4744 


1.331 


4 q94('-3'l 


2 8fi8f-2'l 


-2 184('-3~) 


-1 554('-2~) 


0.80 


4 8n5l'-5'l 


3.684(-4) 


0.3378 


0.8247 


0.6162 


1.505 


8 547('-3'l 


4 3n2('-2'l 


-A 245('-3'l 


-2 64n('-2~) 


0.90 


9.917(-5) 


6.801(-4) 


0.4249 


0.9159 


0.7750 


1.671 


1.374(-2) 


6 16nf-2'l 


-7 624f-3'l 


-4 2n7f-2~) 


1.00 


1.904(-4) 


1.184(-3) 


n 5208 


1.003 


0.9501 


1.829 


2 103^-2") 


8 508('-2"l 


-1 286f-2~) 


-6 377('-2~) 


1.10 


3.453(-4) 


1 968 ('-3'1 


n 6953 


1.085 


1.140 


1.978 




0.1141 


-2 069(^-9'! 


-9 977('-9') 


1.20 


5.973(-4) 


3.151(-3) 


0.7377 


1.162 


1.345 


2.116 


4.407(-2) 


0.1496 


-3 1 7nf-2'l 


-0.1305 


1.30 


9.942(-4) 


4 895 ('-3'1 


0.8575 


1.234 


1.563 


2.243 


6 in9f-2~) 


0.1922 


-4 705^-2") 


-0.1783 


1.40 


1 602 


7.419(-3) 


0.9843 


1.299 


1.793 


2.357 


8 279^-2^ 


0.2432 


-6.774(-2) 


-0.2376 


1.50 


2.514(-3) 


1.102(-2) 


1.117 


1.358 


2.034 


2.456 


0.1100 


3036 


-9 501 f-2^ 


-0.3101 


1.60 


S 856^-^^ 

r_} . \_t I Of 


1.612(-2) 


1.256 


1.411 


2.284 


2.539 


0.1439 


0.3750 


-0.1302 


-0.3970 


1.70 


5 8n6('-3~) 

'J . O yj VJ 1 rj 1 


2 329('-2^ 


1.399 


1.456 


2.542 


2.604 


0.1855 


0.4594 


-0.1749 


-0.4997 


1.80 


8 610f-3~) 


3 335/-2'l 


1.547 


1.494 


2.804 


2.648 


2363 


5593 


-0.2307 


-0.6194 


1.90 


1.261 (-2) 


4.753(-2) 


1.698 


1.526 


3.070 


2.668 


2980 


0.6782 


-0.2994 


-0.7571 


9 nn 


1 QQ1 / 0\ 
± .OOl \~^J 




± .ou^ 


1 "I'll 


o.oo / 






n 891 4 




n Qi "id 

-u.y ±o^ 


2.10 


2.644(-2) 


9.663(-2) 


2.008 


1.572 


3.601 


2.621 


0.4633 


0.9966 


-0.4827 


-1.088 


2.20 


3.809(-2) 


0.1393 


2.166 


1.592 


3.860 


2.544 


0.5735 


1.217 


-0.6010 


-1.281 


2.30 


5.505(-2) 


0.2045 


2.326 


1.620 


4.109 


2.418 


0.7089 


1.505 


-0.7394 


-1.489 


2.40 


8.032(-2) 


0.3100 


2.491 


1.672 


4.342 


2.231 


0.8782 


1.906 


-0.8992 


-1.708 


2.50 


0.1197 


0.4979 


2.663 


1.797 


4.552 


1.956 


1.097 


2.528 


-1.081 


-1.929 


2.60 


0.1867 


0.9026 


2.857 


2.158 


4.728 


1.538 


1.402 


3.713 


-1.284 


-2.130 


2.70 


0.3305 


2.398 


3.133 


3.903 


4.849 


0.7890 


1.917 


7.685 


-1.504 


-2.243 


2.75270 


0.6770 


oo 


3.641 


oo 


4.868 


-0.6130 


2.888 


oo 


-1.620 


-1.969 
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TABLE VI. The functions and their derivatives to represent the equihbrium figure for n = 1. (II) 



n = 1.0 


s 








d'-'^^tb', /d£ 




d''*'i/)99 IdP 




d^^^xjji/di. 




d'^^^xb^, ld£ 


0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


0.000 


0.000 


2 1 7nC-.3"> 


4 336 ('-2'! 


4 281 C-3'1 


8 556 C-2') 


1 193C-6"> 


4.770(-5) 


-5 2n3C-8') 


-2 601 ('-6'1 


0.20 


0.000 


0.000 


8.660(-3) 


8.635(-2) 


1.709(-2) 


0.1704 


1.906(-5) 


3.808(-4) 


-1.663(-6) 


-4.155(-5) 


0.30 


0.000 


0.000 


1.942(-2) 


0.1286 


3 831f-2") 


2538 


9 6271^-51 


1 2811^-31 


-1 260f-5~) 


-2 098f-4'l 


0.40 


0.000 


0.000 


3 434l'-2') 


0.1697 


6 779('-2~) 


3352 


3 033 ('-4^ 


3 022('-3'l 


-5 297('-5~> 


-6 6n5('-4~) 


0.50 


0.000 


0.000 


5.332(-2) 


0.2094 


0.1053 


0.4138 


7.374(-4) 


5.866(-3) 


-1.611(-4) 


-1.605(-3) 


0.60 


0.000 


0.000 


7.617(-2) 


0.2473 


0.1504 


0.4890 


1 521 /-3"1 


1 006/-2"! 


-3 992^^-4^ 


-3 308(^-3^ 


0.70 


0.000 


0.000 


0.1027 


0.2831 


9099 


0.5601 


2.802(-3) 


1 5S3('-9'l 


-8 584('-4') 


-6 n85('-3'; 


0.80 


0.000 


0.000 


0.1327 


0.3165 


9693 


0.6267 


4: . ( '-to \ J 


9 339('-9'l 


-1 664('-3'; 




0.90 


0.000 


0.000 


0.1659 


0.3471 


0.3281 


6880 


7 546f-3'l 


3 292^-2') 


-2 979f-3'l 


-1 634f-2~) 


1.00 


0.000 


0.000 


2n2n 


0.3747 


0.3997 


0.7438 


1.140(-2) 


4 456('-9'l 


-5 nn8('-3') 


-9 465('-9') 


1 in 


n nnn 


n nnn 


n 94ns 


n ^009 


n 47fifi 

U.I 1 uu 


n 7Q^4 






Q nnni' 




1.20 


0.000 


0.000 


0.2817 


0.4201 


5589 


n 8365 


9 31 7(-9) 




_1 99P^('_9^ 
±..ii.iitj\ .ii ( 


-A QQ1 ('-9') 


1.30 


0.000 


0.000 


0.3247 


0.4375 


0.6437 


0.8728 


3 1 55('-9'l 


9 330l'-9'l 


-1 8in('-9') 


-6 780l'-9~) 


1.40 


0.000 


0.000 


0.3691 


0.4511 


0.7325 


0.9019 


4.191(-2) 


0.1143 


-9 5q5('-9') 


-8 Q8fi(^-9'1 


1.50 


0.000 


n nnn 


n 4147 


0.4607 


n 8938 


0.9237 


5 450f-9'l 


0.1378 


-3 6931^-9') 


-0.1165 


1.60 


0.000 


0.000 


0.4611 


0.4664 


0.9170 


n Q3Sn 




0.1636 


-4 q43('-9') 


-0.1483 


1.70 


0.000 


0.000 


0.5079 


0.4681 


1.011 


0.9446 


8 798f'-9'l 


0.1916 


-6 6n7('-9'l 


-0.1856 


1.80 


0.000 


0.000 


0.5546 


0.4656 


1.106 


0.9436 


0.1079 


0.2218 


-8 673('-9~) 


-0.2287 


1.90 


0.000 


0.000 


n fionQ 


0.4591 


1.200 


n Q35n 


0.1317 


n 953Q 


-0.1120 


-0.2781 


2.00 


0.000 


0.000 


0.6463 


0.4485 


1.292 


0.9189 


0.1588 


0.2878 


-0.1426 


-0.3340 


2.10 


0.000 


0.000 


n 6905 


0.4340 


1.383 


n 8955 


0.1893 


n 3933 


-0.1791 


-0.3967 


2.20 


0.000 


0.000 


0.7330 


0.4157 


1.471 


0.8649 


9935 


n 3600 


-0.2221 


-0.4663 


2.30 


n nnn 


n nnn 








n 897fi 


0.2614 




_n 979c; 




2.40 


0.000 


0.000 


0.8116 


0.3681 


1.637 


0.7838 


0.3030 


0.4361 


-0.3310 


-0.6265 


2.50 


0.000 


0.000 


0.8470 


0.3391 


1.713 


0.7341 


0.3486 


0.4748 


-0.3981 


-0.7170 


2.60 


0.000 


0.000 


0.8793 


0.3071 


1.783 


0.6787 


0.3980 


0.5134 


-0.4746 


-0.8142 


2.70 


0.000 


0.000 


0.9083 


0.2722 


1.848 


0.6184 


0.4513 


0.5516 


-0.5611 


-0.9179 


2.80 


0.000 


0.000 


0.9337 


0.2347 


1.907 


0.5537 


0.5083 


0.5889 


-0.6584 


-1.028 


2.90 


0.000 


0.000 


0.9552 


0.1949 


1.959 


0.4851 


0.5690 


0.6250 


-0.7669 


-1.143 


3.00 


0.000 


0.000 


0.9726 


0.1532 


2.004 


0.4135 


0.6332 


0.6594 


-0.8872 


-1.264 


3.10 


0.000 


0.000 


0.9858 
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TABLE VII. The functions and their derivatives to represent the equilibrium figure for n = 1.5. (II) 



n = 1.5 



f 

s 










d^^^^i Idi 








d^'^^ip4/di 






0.00 
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2.10 


/I 
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-2) 


0.2561 


0.1353 


0.5131 


0.3175 


7.448(-2) 


0.1062 
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2.20 


-o 
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0.2691 
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0.5445 


0.3103 
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91 
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0.1147 
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-o 
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0.1090 


0.1220 
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-0.2179 
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-1 
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-2) 


0.3017 
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0.6342 


0.2877 


0.1213 
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2.60 


-1 
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-2) 


-2.646 
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0.3101 
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0.2810 


0.1338 


0.1259 


-0.1649 
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2.70 


-1 
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-3.106 
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0.3171 
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0.6904 


0.2753 


0.1464 


0.1252 


-0.1956 
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2.80 


-1 


932( 


-2) 


-3.634 


-2) 


0.3227 


4 839('-2'l 


0.7177 


0.2710 


0.1588 


0.1224 


-0.2304 


-0.3704 


2.90 


-2 


325( 


-2) 


-4.242 


-2) 


0.3268 


3.291(-2) 


0.7447 


0.2683 


0.1708 


0.1170 


-0.2698 


-0.4181 


3.00 


-2 


784 


-2) 


-4.950 


-2) 


0.3293 


1.656(-2) 


0.7714 


0.2678 


0.1821 


0.1086 


-0.3142 


-0.4703 


3.10 


-3 


319 


-2) 


-5.784 


-2) 


0.3300 


-9.219(-4) 


0.7983 


0.2698 


0.1924 


9.651(-2) 


-0.3641 


-0.5275 


3.20 


-3 
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-2) 


-6.785 


-2) 


0.3290 


-1.994(-2) 
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0.2013 
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-4 
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-8.021 
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0.3260 


-4.126(-2) 


0.8533 


0.2832 
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-5 
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0.3206 


-6.634(-2) 
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3.50 


-6 
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-2) 


-0.1181 


0.3125 


-9.850(-2) 
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0.3140 


0.2124 
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3.60 


-7 
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-2) 


-0.1548 


0.3004 


-0.1493 
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0.3395 


0.2062 
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-0.9174 
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-8 
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-2) 


-0.2311 


0.2908 


-0.2555 


0.9640 


0.3585 
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-0.3006 


-0.7674 


-0.9761 
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TABLE VIII. The functions and their derivatives to represent the equiUbrium figure for n = 2. (II) 



n = 2.0 
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-3.427 


-2) 


0.4182 


0.1982 


6.317 


-2) 


-2.004 


;-2) 


-0.3738 


-0.4438 


4.10 


-4.765( 


-2) 


-4.926 


-2) 


9.080(-2) 


-3.945 


-2) 


0.4388 


0.2127 


6.072 


-2) 


-2.911 


;-2) 


-0.4205 


-0.4900 


4.20 


-5.280( 


-2) 


-5.366( 


-2) 


8.658(-2) 


-4.502 


-2) 


0.4608 


0.2284 


5.731 


-2) 


-3.930 


[-2) 


-0.4720 


-0.5404 


4.30 


-5.840 


-2) 


-5.846 


-2) 


8.178(-2) 


-5.103 


-2) 


0.4845 


0.2456 


5.282 


-2) 


-5.074 


[-2) 


-0.5287 


-0.5954 


4.35287 


-6.156 


-2) 


-6.119 


-2) 


7.899(-2) 


-5.442 


-2) 


0.4977 


0.2552 


4.996 


-2) 


-5.735 


[-2) 


-0.5610 


-0.6264 
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TABLE IX. The velocity potentials and their derivatives for n = 0.5. 



n = 0.5 

















0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


1.646 (-3) 


3.291(-2) 


1.439(-5) 


4.324(-4) 


-1.148(-7) 


-4.623(-6) 


0.20 


6.583(-3) 


6.584(-2) 


1.149(_4) 


1.725(-3) 


-1.828(-6) 


-3.662(-5) 


0.30 


1.481(-2) 


9.878(-2) 


3. 879 (-4) 


3.881(-3) 


-9.248(-6) 


-1.234(-4) 


0.40 


2.634(-2) 


0.1318 


9.196(-4) 


6.901(-3) 


-2.923(-5) 


-2.926(-4) 


0.50 


4.117(-2) 


0.1648 


1.797(-3) 


1.079(-2) 


-7.139(-5) 


-5.717(-4) 


0.60 


5.930(-2) 


0.1978 


3.106(-3) 


1.555(-2) 


-1.481(-4) 


-9.888(-4) 


0.70 


8.074(-2) 


0.2310 


4.936(-3) 


2.119(-2) 


-2.746(-4) 


-1.572(-3) 


0.80 


0.1055 


0.2642 


7.373(-3) 


2.771(-2) 


-4.689(-4) 


-2.350(-3) 


0.90 


0.1336 


0.2975 


1.051(-2) 


3.511(-2) 


-7.519(-4) 


-3.351(-3) 


1.00 


0.1650 


0.3309 


1.442(-2) 


4.341(-2) 


-1.147(-3) 


-4.605(-3) 


1.10 


0.1998 


0.3644 


1.922(-2) 


5.262(-2) 


-1.682(-3) 


-6.141(-3) 


1.20 


0.2379 


0.3981 


2.498(-2) 


6.275(-2) 


-2.386(-3) 


-7.991(-3) 


1.30 


0.2794 


0.4318 


3.180(-2) 


7.379(-2) 


-3.292(-3) 


-1.019(-2) 


1.40 


0.3243 


0.4658 


3.977(-2) 


8.578(-2) 


-4.435(-3) 


-1.276(-2) 


1.50 


0.3725 


0.4999 


4.899(-2) 


9.872(-2) 


-5.856(-3) 


-1.574(-2) 


1.60 


0.4242 


0.5343 


5.955(-2) 


0.1126 


-7.597(-3) 


-1.916(-2) 


1.70 


0.4794 


0.5688 


7.155(-2) 


0.1275 


-9.705(-3) 


-2.307(-2) 


1.80 


0.5380 


0.6036 


8.509(-2) 


0.1435 


-1.223(-2) 


-2.749(-2) 


1.90 


0.6001 


0.6386 


0.1003 


0.1604 


-1.522(-2) 


-3.246(-2) 


2.00 


0.6658 


0.6740 


0.1172 


0.1784 


-1.874(-2) 


-3.804(-2) 


2.10 


0.7349 


0.7097 


0.1360 


0.1975 


-2.285(-2) 


-4.425(-2) 


2.20 


0.8077 


0.7458 


0.1567 


0.2178 


-2.761(-2) 


-5.115(-2) 


2.30 


0.8841 


0.7823 


0.1796 


0.2392 


-3.310(-2) 


-5.880(-2) 


2.40 


0.9642 


0.8192 


0.2046 


0.2618 


-3.940(-2) 


-6.723(-2) 


2.50 


1.048 


0.8568 


0.2320 


0.2858 


-4.658(-2) 


-7.653(-2) 


2.60 


1.136 


0.8949 


0.2618 


0.3111 


-5.473(-2) 


-8.676(-2) 


2.70 


1.227 


0.9339 


0.2943 


0.3380 


-6.396(-2) 


-9.802(-2) 


2.75270 


1.277 


0.9548 


0.3125 


0.3528 


-6.930(-2) 


-0.1044 
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TABLE X. The velocity potentials and their derivatives for n = 1. 



n = 1.0 


s 














0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


9 819f-4') 


1 964('-2') 


7 695(^-6") 


2.313(-4) 


-5 453C-8") 


-2 195('-6') 


0.20 


3.929(-3) 


3.930(-2) 


6.150(-5) 


9.233(-4) 


-8.684(-7) 


-1.740(-5) 


0.30 


8.844(-3) 


5.901(-2) 


2.077(-4) 


2.079(-3) 


-4.396(-6) 


-5.871(-5) 


0.40 


1 573(-2) 


7 878f-2) 


4.927(-4) 


3 702f-3l 


-1.391(-5) 


-1.394(-4) 


0.50 


2.460(-2) 


9.863(-2) 


9.636(-4) 


5.797(-3) 


-3.401(-5) 


-2.728(-4) 


0.60 


3.546(-2) 


0.1186 


1.668(-3) 


8.371(-3) 


-7.068(-5) 


-4.729(-4) 


0.70 


4 832(-2) 


0.1387 


2 654f-3l 


1.143(-2) 


-1.313(-4) 


-7.539(-4) 


0.80 


6 320f-2~) 


0.1589 


3.971(-3) 


1.499(-2) 


-2.246(-4) 


-1.130(-3) 


0.90 


8 012f-2~) 


0.1794 


5 668f-3l 

^ ^yjyjy_y \ 


1 905 ('-21 


-3.610(-4) 


-1.618(-3) 


1.00 


9 908f-2~) 


2000 


7.799(-3) 


2.364(-2) 


-5 523f-4l 


-2.232(-3) 


1.10 


0.1201 


0.2209 


1.041(-2) 


2 876f-2l 


-8.121(-4) 


-2.989(-3) 


1.20 


0.1433 


0.2420 


1.357(-2) 


3.444(-2) 


-1.155(-3) 


-3 q08f-3~) 


1.30 


0.1685 


0.2634 


1.732(-2) 


4 0691^-21 


-1 600('-3l 


-5 007f-3~) 


1.40 


0.1960 


0.2851 


2.173(-2) 


4.753(-2) 


-2.164(-3) 


-6 306f-3~) 


1.50 


0.2256 


0.3072 


2.685(-2) 


5.499(-2) 


-2.868(-3) 


-7.827(-3) 


1.60 


0.2574 


3296 


3 275 f-2") 


6 310f-2l 


-3.737(-3) 


-9 5Ql('-3~) 


1.70 


0.2915 


0.3525 


3.949(-2) 


7.188(-2) 


-4.795(-3) 


-1.162(-2) 


1.80 


0.3279 


0.3758 


4.714(-2) 


8.136(-2) 


-6.072(-3) 


-1.395(-2) 


1.90 


0.3667 


3996 


5 579(^-2^ 


9 159(^-2^ 


-7.597(-3) 


-1.661(-2) 


2.00 


0.4079 


0.4240 


6 549f-2l 


0.1026 


-9 405(^-3^ 


-1.961(-2) 


2.10 


0.4515 


0.4490 


7 633f-2l 


0.1145 


-1.153(-2) 


-2.301(-2) 


2.20 


0.4977 


0.4745 


8.841(-2) 


0.1272 


-1.402(-2) 


-2.684(-2) 


2.30 


0.5464 


5008 


0.1018 


0.1408 


-1.692(-2) 


-3.113(-2) 


2.40 


0.5979 


0.5279 


0.1166 


0.1555 


-2.026(-2) 


-3.593(-2) 


2.50 


0.6520 


0.5558 


0.1329 


0.1712 


-2.412(-2) 


-4.130(-2) 


2.60 


0.7090 


0.5845 


0.1509 


0.1880 


-2.854(-2) 


-4.728(-2) 


2.70 


0.7690 


0.6143 


0.1706 


0.2061 


-3.360(-2) 


-5.394(-2) 


2.80 


0.8319 


0.6451 


0.1922 


0.2255 


-3.936(-2) 


-6.135(-2) 


2.90 


0.8980 


0.6771 


0.2157 


0.2463 


-4.590(-2) 


-6.957(-2) 


3.00 


0.9674 


0.7104 


0.2415 


0.2687 


-5.330(-2) 


-7.871(-2) 


3.10 


1.040 


0.7450 


0.2695 


0.2927 


-6.167(-2) 


-8.884(-2) 


3.111")!) 


1.071 


().7r)!)!) 


().2cSl<) 


().:i():!2 


-G.51()(-2) 


-!).:!:!7(-2) 
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TABLE XI. The velocity potentials and their derivatives for n = 1.5. 



n = 1.5 


s 














0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


5.484(-4) 


1.097(-2) 


3.754(-6) 


1.129(-4) 


-2.308(-8) 


-9.293(-7) 


0.20 


2.195(-3) 


2.196(-2) 


3.001(-5) 


4.507(-4) 


-3.677(-7) 


-7.372(-6) 


0.30 


4.943(-3) 


3.301(-2) 


1.014(-4) 


1.016(-3) 


-1.863(-6) 


-2.490(-5) 


0.40 


8.798(-3) 


4.412(-2) 


2.409(-4) 


1.812(-3) 


-5.902(-6) 


-5.923(-5) 


0.50 


1.377(-2) 


5.533(-2) 


4.717(-4) 


2.844(-3) 


-1.446(-5) 


-1.162(-4) 


0.60 


1.987(-2) 


6.666(-2) 


8.177(-4) 


4.118(-3) 


-3.010(-5) 


-2.021(-4) 


0.70 


2.711(-2) 


7 813f-2~) 


1 303f-3l 


5.641(-3) 


-5 603f-5^ 


-3.232(-4) 


0.80 


3 550f-2~) 


8 978f-2) 


1.954(-3) 


7.422(-3) 


-9 613f-5l 


-4.865(-4) 


0.90 


4 507('-2~) 


0.1016 


2.797(-3) 


9.473(-3) 


-1.550(-4) 


-6.993(-4) 


1.00 


5 583f-2~) 


0.1137 


3 858f-3l 


1.181(-2) 


-2.379(-4) 


-9.695(-4) 


1.10 


6 782f-2~) 


0.1260 


5 168(^-31 


1.443(-2) 


-3.510(-4) 


-1.306(-3) 


1.20 


8.105(-2) 


0.1386 


6.755(-3) 


1.737(-2) 


-5.015(-4) 


-1.717(-3) 


1.30 


9 555f-2~) 


0.1515 


8 6531^-31 


2.064(-2) 


-6.972(-4) 


-2.213(-3) 


1.40 


0.1114 


0.1648 


1 089f-2"l 


2.425(-2) 


-9.473(-4) 


-2 806('-3~) 


1.50 


0.1285 


0.1784 


1.352(-2) 


2.824(-2) 


-1.262(-3) 


-3 506f-3~) 


1.60 


0.1471 


0.1925 


1 655 f-2") 


3.262(-2) 


-1 653f-3l 


-4.328(-3) 


1.70 


0.1670 


0.2070 


2.005(-2) 


3.741(-2) 


-2.132(-3) 


-5.285(-3) 


1.80 


0.1885 


0.2219 


2.405(-2) 


4.266(-2) 


-2.715(-3) 


-6.394(-3) 


1.90 


0.2114 


0.2374 


2 860(^-2^ 


4 838(^-2^ 


-3.416(-3) 


-7.671(-3) 


2.00 


n 2360 


n 2535 


3 375 ('-2") 


5 462^-2") 


-4 255f-3") 


-9.136(-3) 


2.10 


0.2621 


0.2701 


3 954f-2l 


6.140(-2) 


-5 250f-3l 


-1.081(-2) 


2.20 


2900 


0.2874 


4 605f-2l 


6 878f-2l 


-6 425f-3l 


-1.271(-2) 


2.30 


0.3196 


0.3054 


5 332f-2l 


7.679(-2) 


-7 802f-3l 


-1.487(-2) 


2.40 


0.3511 


0.3241 


6.143(-2) 


8.547(-2) 


-9 409(^-3^ 


-1.732(-2) 


2.50 


0.3845 


0.3436 


7.044(-2) 


9 488f-2l 


-1.128(-2) 


-2.008(-2) 


2.60 


0.4199 


3639 


8 043(^-21 


0.1051 


-1.344(-2) 


-2.318(-2) 


2.70 


0.4573 


0.3851 


9.148(-2) 


0.1161 


-1.592(-2) 


-2.666(-2) 


2.80 


0.4969 


0.4073 


0.1037 


0.1280 


-1.878(-2) 


-3.057(-2) 


2.90 


0.5388 


0.4305 


0.1171 


0.1409 


-2.205(-2) 


-3.493(-2) 


3.00 


0.5830 


0.4547 


0.1319 


0.1548 


-2.579(-2) 


-3.981(-2) 


3.10 


0.6298 


0.4801 


0.1481 


0.1698 


-3.003(-2) 


-4.523(-2) 


3.20 


0.6791 


0.5066 


0.1659 


0.1860 


-3.485(-2) 


-5.127(-2) 


3.30 


0.7311 


0.5343 


0.1854 


0.2035 


-4.031(-2) 


-5.797(-2) 


3.40 


0.7860 


0.5634 


0.2066 


0.2223 


-4.647(-2) 


-6.539(-2) 


3.50 


0.8438 


0.5938 


0.2299 


0.2424 


-5.341(-2) 


-7.360(-2) 


3.60 


0.9048 


0.6256 


0.2552 


0.2641 


-6.122(-2) 


-8.266(-2) 


3.65375 


0.9389 


0.6433 


0.2697 


0.2764 


-6.580(-2) 


-8.790(-2) 



51 



TABLE XII. The velocity potentials and their derivatives for n = 2. 



n = 2.0 


s 














0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


2.801 (-4) 


5 6QA(-'S) 


1.632(-6) 


4.909(-5) 


-8.494(-9) 


-3.420(-7) 


0.20 


1.121(-3) 


1.123(-2) 


1.306(-5) 


1.962 (-4) 


-1.354(-7) 


-2.716(-6) 


0.30 


2 527f-3) 


1 689f-2^ 


4.416(-5) 


4.430(-4) 


-6.869(-7) 


-9.190(-6) 


0.40 


4 502(-3) 


2.261(-2) 


1.050(-4) 


7.917(-4) 


-2.179(-6) 


-2.190(-5) 


0.50 


7.052(-3) 


2.842(-2) 


2.060(-4) 


1.246(-3) 


-5.347(-6) 


-4.310(-5) 


0.60 


1.019(-2) 


3.432(-2) 


3.577(-4) 


1.809(-3) 


-1.116(-5) 


-7.519(-5) 


0.70 


1.392(-2) 


4 034f-2~) 


5.715(-4) 


2.487(-3) 


-2 083f-5l 


-1.207(-4) 


0.80 


1 S26('-2~) 


4 65n('-2~) 


8.591(-4) 


3 2Sfi('-3'l 


-3.584(-5) 


-1.825(-4) 


0.90 


2.323(-2) 


5 283f-2~) 


1.233(-3) 


4.212(-3) 


-5.797(-5) 


-2.637(-4) 


1.00 


2 883f-2~) 


5 933f-2'l 


1.706(-3) 


5.276(-3) 


-8 933f-5') 


-3.675(-4) 


1.10 


3 510f-2~) 


6 605f-2') 


2.293(-3) 


6.485(-3) 


-1.324(-4) 


-4.978(-4) 


1.20 


4.205(-2) 


7 300f-2~) 


3 008f-3') 


7.850(-3) 


-1.899(-4) 


-6.587(-4) 


1.30 


4.971(-2) 


8 020('-2~l 


3 868(^-31 


9.384(-3) 


-2.653(-4) 


-8.548(-4) 


1.40 


5 810('-2i 


8 767(^-2) 


4.891 (-3) 


1.110(-2) 


-3.622(-4) 


-1.091(-3) 


1.50 


6 725f-2^ 


9 545f-2l 


6 095f-3l 


1.301(-2) 


-4.850(-4) 


-1.373(-3) 


1.60 


7.720(-2) 


0.1036 


7.500(-3) 


1.513(-2) 


-6.386(-4) 


-1.708(-3) 


1.70 


8.798(-2) 


0.1120 


9.128(-3) 


1.748(-2) 


-8.285(-4) 


-2.101(-3) 


1.80 


9 962('-2~) 


0.1209 


1.100(-2) 


2.008(-2) 


-1.061 (-3) 


-2.561(-3) 


1.90 


0.1122 


0.1301 


1.315(-2) 


2.294(-2) 


-1.343(-3) 


-3.097(-3) 


2.00 


0.1257 


0.1398 


1.560(-2) 


2 609^-2") 


-1.683(-3) 


-3.717(-3) 


2.10 


0.1401 


0.1500 


1.838(-2) 


2.955(-2) 


-2 090f-3l 


-4.432(-3) 


2.20 


0.1557 


0.1607 


2.152(-2) 


3 3351^-21 


-2.573(-3) 


-5.252(-3) 


2.30 


0.1723 


0.1718 


2.506(-2) 


3.751(-2) 


-3.144(-3) 


-6.192(-3) 


2.40 


0.1901 


0.1836 


2.904(-2) 


4.206(-2) 


-3.816(-3) 


-7.263(-3) 


2.50 


2090 


0.1960 


3.349(-2) 


4.703(-2) 


-4.602(-3) 


-8.480(-3) 


2.60 


0.2293 


2089 


3.846(-2) 


5.245(-2) 


-5.517(-3) 


-9 860f-3") 


2.70 


2508 


0.2226 


4.399(-2) 


5 835f-2l 


-6.579(-3) 


-1.142(-2) 


2.80 


0.2738 


2369 


5.014(-2) 


6.477(-2) 


-7.808(-3) 


-1.318(-2) 


2.90 


0.2983 


0.2520 


5.697(-2) 


7.175(-2) 


-9.222(-3) 


-1.515(-2) 


3.00 


0.3242 


0.2678 


6.451(-2) 


7.932(-2) 


-1.085(-2) 


-1.737(-2) 


3.10 


0.3518 


0.2844 


7.285(-2) 


8.753(-2) 


-1.270(-2) 


-1.985(-2) 


3.20 


0.3811 


0.3019 


8.204(-2) 


9.641(-2) 


-1.483(-2) 


-2.261(-2) 


3.30 


0.4122 


0.3202 


9.216(-2) 


0.1060 


-1.724(-2) 


-2.570(-2) 


3.40 


0.4452 


3395 


0.1033 


0.1164 


-1.998(-2) 


-2.912(-2) 


3.50 


0.4802 


3596 


0.1155 


0.1275 


-2.307(-2) 


-3.291(-2) 


3.60 


0.5172 


0.3808 


0.1288 


0.1395 


-2.657(-2) 


-3.711(-2) 


3.70 


0.5564 


0.4029 


0.1434 


0.1524 


-3.051(-2) 


-4.174(-2) 


3.80 


0.5978 


0.4261 


0.1593 


0.1663 


-3.494(-2) 


-4.685(-2) 


3.90 


0.6416 


0.4503 


0.1767 


0.1812 


-3.990(-2) 


-5.246(-2) 


4.00 


0.6879 


0.4756 


0.1956 


0.1971 


-4.544(-2) 


-5.861(-2) 


4.10 


0.7368 


0.5020 


0.2161 


0.2140 


-5.164(-2) 


-6.535(-2) 


4.20 


0.7883 


0.5296 


0.2384 


0.2322 


-5.853(-2) 


-7.270(-2) 


4.30 


0.8427 


0.5583 


0.2626 


0.2515 


-6.620(-2) 


-8.073(-2) 


4.35287 


0.8726 


0.5739 


0.2762 


0.2622 


-7.059(-2) 


-8.525(-2) 
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TABLE XIII. The normalized radii for spherical stars, the orbital separations at the contact point, the energies, the reduced 
quadrupole moments, and the changes in the central density for n = 0.5, 1, 1.5 and 2. 



n 




Rc/o-o 


Eself 






*ii/(Ma§) 


5pc/Pc 


0.5 


2.75270 


2.496 


-1.111 


-0.5000 


1.962 


0.6691 


2.069 


1.0 


3.14159 


2.426 


-1.000 


-0.5000 


1.040 


0.3465 


2.280 


1.5 


3.65375 


2.384 


-0.8571 


-0.5000 


0.5749 


0.1915 


2.158 


2.0 


4.35287 


2.358 


-0.6667 


-0.5000 


0.2958 


9.861 (-2) 


1.924 



TABLE XIV. The individual energies for n = 0.5, 1, 1.5 and 2. 



n 


Ilsei/ 






{Wself)quad 


i^^A^int) point 


(^'^irit ) quad 


Tpoint 


Tquad 


0.5 


0.2222 


-0.9946 


-1.333 


1.953 


-1.000 


-2.007 


0.5000 


3.011 


1.0 


0.5000 


-1.300 


-1.500 


1.819 


-1.000 


-1.040 


0.5000 


1.559 


1.5 


0.8571 


-1.437 


-1.714 


1.724 


-1.000 


-0.5746 


0.5000 


0.8620 


2.0 


1.333 


-1.479 


-2.000 


1.627 


-1.000 


-0.2958 


0.5000 


0.4437 
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TABLE XV. The total energy, the total angular momentum and the orbital angular velocity along the equilibrium sequences 
of the binary system. 



R/ao 


E/{M^/ao) 




J/{M^aoY'^ 


n/(7rpo)^/' 


n = 0.5 


7.0 


-1.183 




1.871 


8.819(-2) 


6.5 


-1.188 




1.803 


9.857(-2) 


6.0 


-1.194 




1.733 


0.1112 


5.5 


-1.202 




1.659 


0.1267 


5.0 


-1.211 




1.583 


0.1462 


4.5 


-1.222 




1.502 


0.1713 


4.0 


-1.236 




1.418 


0.2047 


3.5 


-1.253 




1.330 


0.2508 


3.0 


-1.275 




1.240 


0.3181 


2.5 


-1.303 




1.152 


0.4257 


2.496 


-1.303 




1.151 


0.4269 


n = 1 


7.0 


-1.071 




1.871 


8.818(-2) 


6.5 


-1.077 




1.803 


9.855(-2) 


6.0 


-1.083 




1.732 


0.1111 


5.5 


-1.091 




1.659 


0.1266 


5.0 


-1.100 




1.582 


0.1461 


4.5 


-1.111 




1.501 


0.1712 


4.0 


-1.125 




1.416 


0.2044 


3.5 


-1.142 




1.327 


0.2501 


3.0 


-1.165 




1.233 


0.3163 


2.5 


-1.196 




1.136 


0.4197 


2.426 


-1.201 




1.122 


0.4400 


n = 1.5 


7.0 


-0.9286 




1.871 


8.818(-2) 


6.5 


-0.9341 




1.803 


9.855(-2) 


6.0 


-0.9405 




1.732 


0.1111 


5.5 


-0.9480 




1.659 


0.1266 


5.0 


-0.9571 




1.582 


0.1461 


4.5 


-0.9682 




1.501 


0.1711 


4.0 


-0.9820 




1.415 


0.2043 


3.5 


-0.9997 




1.325 


0.2498 


3.0 


-1.023 




1.229 


0.3154 


2.5 


-1.055 




1.128 


0.4167 


2.384 


-1.064 




1.104 


0.4487 


n = 2 


7.0 


-0.7381 




1.871 


8.818(-2) 


6.5 


-0.7436 




1.803 


9.854(-2) 


6.0 


-0.7500 




1.732 


0.1111 


5.5 


-0.7576 




1.658 


0.1266 


5.0 


-0.7666 




1.581 


0.1461 


4.5 


-0.7777 




1.500 


0.1711 


4.0 


-0.7916 




1.415 


0.2042 


3.5 


-0.8094 




1.324 


0.2496 


3.0 


-0.8329 




1.227 


0.3148 


2.5 


-0.8654 




1.123 


0.4150 


2.358 


-0.8770 




1.092 


0.4539 
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FIG. 1. (a) The total energy and (b) the total angular momentum as functions of the orbital separation, and (c) the orbital 
angular velocity as functions of the total angular momentum for the polytropic index n — 0.5. Solid line denotes our results. 
Filled circles and open triangles are the results of Lai, Rasio and Shapiro (1994) and Uryu and Eriguchi (1998). 
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FIG. 2. (a) The total energy and (b) the total angular momentum as functions of the orbital separation, and (c) the orbital 
angular velocity as functions of the total angular momentum for the polytropic index n = 1. Solid line denotes our results. 
Filled circles and open triangles are the results of Lai, Rasio and Shapiro (1994) and Uryu and Eriguchi (1998). 
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FIG. 3. (a) The total energy and (b) the total angular momentum as functions of the orbital separation, and (c) the orbital 
angular velocity as functions of the total angular momentum for the polytropic index n = 1.5. Solid line denotes our results. 
Filled circles and open triangles axe the results of Lai, Rasio and Shapiro (1994) and Uryu and Eriguchi (1998). 
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FIG. 4. (a) The total energy and (b) the total angular momentum as functions of the orbital separation, and (c) the orbital 
angular velocity as functions of the total angular momentum for the polytropic index n = 2. Solid line denotes our results. 
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FIG. 5. The relative errors for the functions presented the deformations of the figure for the polytropic index n = 1. Thick 
soUd, dotted, dashed, long-dashed, dot-dashed and thin solid lines denote the relative errors for ^^Vz, ^^Vsj ''^V4i ''^Vz, '^V4 
and '^Vsi respectively. 
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